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Abstract 

We discuss supersymmetric Yang-Mills theories with the multiple scales in the brane 
language. The issue concerns N — 2 SUSY gauge theories with massive fundamental mat- 
ter including the UV finite case of riy — 2nc, theories involving products of SU{n) gauge 
groups with bifundamental matter, and the systems with several parameters similar to 
^QCD- We argue that the proper integrable systems are, accordingly, twisted XXX SL{2) 
spin chain, SL{p) magnets and degenerations of the spin Calogero system. The issue of 
symmetries underlying integrable systems is addressed. Relations with the monopole sys- 
tems are specially discussed. Brane pictures behind all these integrable structures in the 
IIB and M theory are suggested. We argue that degrees of freedom in integrable systems 
are related to KK excitations in M theory or D-particles in the IIA string theory, which 
substitute the infinite number of instantons in the field theory. This implies the presence 
of more BPS states in the low-energy sector. 
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1 Introduction 



During the last years it has become clear that the low energy sector of A'^ = 2 supersymmetric gauge theories 
jl], ^ is governed by integrable systems with the finite number of degrees of freedom so that physical quantities 
in the gauge theories (couplings and masses) turn out to be parameters of these finite-dimensional systems. 
There have been considered several examples that include the pure gauge theory |^ as well as the theories 
with adjoint |^ and fundamental matter In particular, in |^ the N — 2 super-QCD with the number of 
fiavors n/ < 2n and the gauge group SU{n) |^] has been considered as a first example of the theory with 
several scales treated in the integrability framework. In the present paper, we discuss more examples of the 
asymptotically free theories with the multiple scales. This condition implies the set of possibilities. Indeed, 
one can consider the theory with massive fundamental matter, the theory with many Aqcd type parameters 
instead of the single one, and the system with the gauge group being product of simple factors Q so that 
there is a parameter attached to each factor. We find integrable counterparts to all these examples. 

In fact, there are some (apparently different) models among field theories, supersymmetric monopoles, 
string theory compactifications which possess the same integrable structure that naturally adds to this list 
another integrable system counterpart. The whole picture is drawn in Fig.l. 
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Figure 1: Various aspects (columns) of the same models (rows) 

Let us stress the special role of quiver diagrams in all these models. The same diagram defines also the 
corresponding brane configurations and integrable systems. It also labels the type of singularity in Calabi-Yau 
manifold in the geometric engineering approach as well as S-duality group. 

It is clear that the revealing of the hidden integrable structure in the vacuum sector of SYM theory can 
not be the aim by itself, thus, we have to explain the whole motivation of the approach. First, let us note 
that the appearance of the integrable system can be expected from the very beginning. Indeed, integrable 
system typically has its configuration or phase spaces realized as some moduli space. That is what we have in 
the super- Yang-Mills (SYM) case: in the pure N = 2 gauge theory, the Coulomb moduli space and instanton 
moduli space are involved. Therefore, we would actually expect a pair of integrable systems. This is, indeed, 
the case, and the two integrable systems - the Whitham and Hitchin-like ones are exhibited. Identification of 
the Whitham dynamics as the one on the Coulomb moduli is quite transparent but the direct relation of the 
Hitchin dynamics with the instanton moduli space is still to be clarified. It has been shown that it is related 
to the monopole rather than to instanton dynamics. Still it is easy to connect these two - since the monopole 
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solutions can be constructed via an infinite dimensional counterpart of the instantonic ADHM construction Q, 
one could naturally identify the instantons taken in the conformal t'Hooft anzatz and organized into ric infinite 
chains with the circle of ric monopoles In what follows we use a kind of Born-Oppenheimer approximation 
and freeze the "slow" Whitham dynamics considering only the "fast" Hitchin one. 

The monopoles which is one of the main subjects of the papeij^ have a specific nature - their role is to 
collect all the instanton nonperturbative contributions in 4d theory. That is why the brane configurations we 
discuss necessarily include the objects which provide the nonperturbative (instanton) effects in the field theory 
so that, as a result, DO-branes within D4-brane gives rise to the nonperturbative solutions of the monopole 
type. Note that we add this ingredient to the M theory approach in comparison with Q, which corresponds to 
additional KK degrees of freedom on the M5 worldvolume. This is in agreement with the general description 
of the finite-gap solutions of integrable systems which can be formulated as dynamics of points on the spectral 
curve. Within this approach the Lax equation is responsible for the structure of the fermionic zero modes in the 
monopole background, while the Baker- Akhiezer function would be associated with the properly normalized 
fermionic zero mode itself. The role of the spectral curve is to capture the band structure of the zero modes 
which are delocalized in the topologically nontrivial background. Note that just analogous band structure of 
the zero modes is expected in the usual QCD. Let us remark that there are similarities with the scenario of 
"confining strings" which also are the effective objects which substitute the infinite number of instantons. 

Turn now to the role which integrable systems play in the context of the low-energy SYM theories. First 
of all, the interpretation of the field-theory vacuum dynamics in terms of finite number degrees of freedom 
immediately raises the question, what these degrees of freedom really are. Now they are recognized as reflecting 
the dynamics of the flnite number of interacting branes. Second, they provide the powerful computational tools 
to deal with the vacuum moduli spaces. Although the D-brane approach looks more transparent, it actually 
reflects the degrees of freedom of integrable systems and, therefore, the quantitative analysis available in the 
dynamical system seems to be very useful. For instance, the phase spaces of integrable systems always have 
an algebraic nature, thus, one usually has a huge underlying symmetry group for the corresponding dynamics 
. Typically, they are of Yangian or quantum affine type. It is natural to look for the counterparts of these 
symmetries in the brane language. These algebras actually gives a field theory analog of the string "BPS 
like" algebras. We comment this subject in section 8. Let us remark that, amazingly, this point of view 
sometimes can provide even with a more convenient physical insight . Another important point about the 
integrable treatment is that the finite dimensional integrable systems admit the Lagrangian description that 
can be considered as the starting point for the quantization of the D-brane transverse degrees of freedom. 

Let us note that in this paper we intensively use the term "multiscale" that we attribute to the three 
different cases. First, we consider the N — 2 theory with fundamental matter where scales are introduced by 
the quark masses. The second example deals with the "group product" case with the distance between NS5 
branes taken as dimensional parameters^. And, finally, the set of parameters can enter the low energy sector 
via the multiple regulator parameters due to a proper dimensional transmutation procedure. 

The paper is organized as follows. At first, we review the picture behind the theory with a single scale 
and discuss the relevance of the brane degrees of freedom for integrability. This is done in section 2. In section 
3 we complete the identification of the integrable structure for the SU (n) theory with fundamentals and show 
that the theory with n/ — 2n is governed by the twisted XXX SL{2) spin chain. Section 4 is devoted to 
the clarification of the role of higher magnets as the integrable counterparts for the "group product" family 
of theories. In the course of our investigation, we works mainly with the "lattice" Lax representations 
that can be characterized by the quadratic Poisson brackets with numerical r-matrix. Still, to deal with the 



^Fix right here that, through the paper, speaking of monopoles, we often imply just a solution to the 
Nahm equations without specializing its boundary conditions. These Nahm equations is the result of applying 
the ADHM construction, and, to describe the true monopole, they are to be added by the proper boundary 
conditions. We return to the discussion of this point in s.6. 

^ These parameters can be associated with the finite cut-off in the UV finite theories or with the dimensional 
parameters of AQCD-type in the asymptotically free ones. 
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multiscale theories, we need "large" Lax representation of the spin Calogero type that is characterized by 
the linear Poisson brackets - maybe with dynamical r-matrix (this is, of course, nothing but Hitchin system 
type representation 0, ^). Quite astonishing, in many integrable systems one can meet these both 
Lax representations so that they lead to the same spectral curve. This fact, along with the different Lax 
representations and their interpretation in brane terms is contained in section 5. The fact that higher magnets 
emerge not accidentally gets more clear in section 6, where the correspondence between spin chains and 
monopoles is established. A new type of gauge theories with multiple A type parameters which is described 
by degenerations of the spin Calogero (Toda) system is introduced in Section 7. We also discuss there the 
corresponding brane interpretation. Some comments on the hidden symmetry groups are presented in section 
8. In the Appendix, we demonstrate some very explicit formulas for the SL{3) magnet that are aimed to 
illustrate some statements of sections 4 and 5 concerning the "group product" case. 



2 Single scale theory 
2.1 Generalities 

We start with the general description of the single scale SUSY gauge theories. Hereafter SUSY YM theories 
are treated as the theories on the D-brane world- volumes. The starting point is the consideration of n parallel 
D-branes embedded in higher dimensions. In order to produce 4d theory, this world has 3+1 noncompact 
dimensions while the coordinates of D-branes in other dimensions are identified with the vacuum expectation 
values of the scalar components in the adjoint hypermultiplets of the theory on the world volume. In fact, 
some of these D-brane coordinates in the transverse directions provide the phase space for the integrable 
many-body problems. 

On the field theory side we consider the softly broken A*' = 4 SYM theory. This theory in the N = 1 
language has three adjoint hypermultiplets whose complex scalar components parametrize the remaining six 
dimensions. The potential energy of these scalars contains the term V = ^Tr^^^^ 2 3^''-' 't't]^ coming from 
the kinetic term and the contribution from the superpotential W — Tr<l>i[$2, ^3] + M^Tr("I>2 + ^s)^ where 
$i are the A' = 1 adjoint superfields and M is the mass of the hypermultiplet. Minimization of the potential 
energy provides classical vacuum configurations that, in our case, correspond to the three scalar fields satisfying 
the nontrivial commutation relations. Note that the large mass term in the Lagrangian freezes the fiuctuations 
transverse to the surface 

Tr{<S>l + ^l) ^ const (2.1) 

which is an algebraic curve in the two-dimensional complex space. 

In this section we will try to fix the degrees of freedom of the integrable systems in brane terms. Actually, 
we need to identify coordinates, momenta and coupling constants in integrable many-body systems. We will 
show that the momenta in the integrable system are related to the positions of D4-branes in the proper 
direction. These branes provide the worldvolume for the d=4 N — 2 theory. Coordinates in the integrable 
system come from positions of other branes responsible for the nonperturbative effects [|l^ , while the coupling 
constants correspond to the regulator masses (say, the mass of adjoint scalar in the Calogero case). 

One more ingredient of the integrability to be explained concerns the meaning of the Lax equation and 
the Baker- Akhiezer function. Let us turn, for the sake of simplicity, to the Toda chain case. There are some 
arguments in favor of interpretation of the Baker-Akhiezer function as the fermion zero mode in a monopole 
background. Indeed, the eigenvalue problem acquires the form of the equation for the fermionic zero mode 
"Hamiltonian" and similar to the Peierls model treatment 

Ci$i+i 4- Ci_i$i_i -f pi$i = A$i (2.2) 

where d = exp(gi+i —qt) and boundary condition ^i+„(X) — exp(in7r(A))$i(A) are imposed, 7r(A) being quasi- 
momentum. The elements of the Lax operator admit the interpretation as matrix elements of the spectral 
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parameter operator 



(2.3) 



The spectral curve with this interpretation plays the role of the dispersion law for the delocalized fermion zero 
modes. 

There are two equivalent ways of constructing Seiberg-Witten solutions from string theory. The relation 
between them was discussed in ^] . The first one, called geometric engineering, deals with the Calabi- 

Yau threefold compactification in the vicinity of singularity |l^ . The other one is the field theory on the 
worldvolume of D4-branes stretched between parallel NS5-branes |h]: 




(2.4) 



0123456789 
+ + + + + + -- -- 
+ + + + -- + -- - 
+ + + + -- - + + + 
where plus stands for extended direction of each brane and minus denotes transverse direction. Since fourbranes 
are finite in the direction, the long-range theory is a four-dimensional N = 2 gauge theory as we need. 
In the IIB case, the original supersymmetry (32 real supercharges) is broken down to A'' = 2 in d = 4 (8 
real supercharges) via geometry of the Calabi-Yau manifold. D-branes, wrapped around shrinking cycles, 
correspond to BPS states that further break half of the supersymmetries. In the IIA case, the SUSY is broken 
by inserting branes into the (perturbatively) flat space^. The BPS states in this picture are membranes of 
M theory ending on a single NS5-brane. For the sake of unity that we are trying to reach in this paper, we 
consider both the Type IIA/M theory and the Type IIB theory in parallel. 



2.2 Type IIB perspective 

If one deals with the Type IIB/F theory, 4ci SU{n) theories without matter emerge as the world volume 
theories of the D7-branes wrapped around an elliptically fibered K3. In fact, the whole picture in this case 
is not so clear. Still some of the important ingredients can be observed. In particular, the key ingredient 
of the integrable system - Calogero model (see, e.g., and notations therein) - namely, the bare spectral 
curve (torus) whose modulus is the bare coupling constant enters the picture automatically. Another crucial 
ingredient - holomorphic 1-form "I> that, after all, turns out to be the Lax operator corresponds to the D7- 
branes while a particular evaluation representation at the marked point is represented by the background 
branes. The n D7-branes are wrapped around the bare torus, therefore, their coordinates in the transverse 
directions are captured by the matrix "1>, if the torus is embedded into the two-dimensional complex surface 
as the supersymmetric cycle [^. The degrees of freedom coming from the antiholomorphic form A (whose 
diagonal elements are the coordinates of the Calogero particles) can be attributed to the branes representing 
nonperturbative efl^ects. Finally, the regulators enter the problem via the insertion of a Wilson lineQ in the 
proper representation at the fiber torus in K3. Eigenvalues of the Wilson line provide the scales in the theory. 
Proper treatment of the regulators is still incomplete but some identification in the orientifold terms seems 
attractive. 

Let us emphasize that the role of the 11-12 dimensions in F theory in this consideration is only to provide 
the necessary input of the symmetry structure, that is, affine algebra and its representation as well as the bare 
coupling constant of d=4 theory. Note also that, in this picture, the branes are wrapped around the bare not 
the actual spectral curve. 

^For models without D6-branes. 

"^This Wilson line provides the spin variables in the spin Calogero model [ p^ p^ ] and reduces to just a 
couphng constant in the standard (spinless) Calogero model. 
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To provide some additional feeling about degrees of freedom of the integrable system let us consider the 
almost trivial "free" many-body integrable system. It corresponds to the d=5 theory without adjoint matter 
. The convenient look at this system follows from the Chern-Simons (CS) interpretation of the Ruijsenaars 
system, now without interaction. The proper picture can be presented in terms of the SU(n) bundle over 
torus if one considers the CS theory on the product of this torus and an interval. The relation for monodromies 
around the cycles gA,gB reads 

gAgBQA^QB^ = 1 

and the diagonal elements of matrices gA,gB which, in the free case, can be diagonalized simultaneously 
represent the coordinates and the momenta of the integrable system gA = diag{e"'); gB = diag{e^''''). Remark 
that, from the F theory point of view, we have degenerated the fiber torus to a cylinder. 

Let us compare this data with the brane picture. There are two coordinates in the brane picture relevant 
for the d=5 theory (let us denote them x"',x^'), the typical size of the brane configuration in the direction 
being R^^ and that in the direction - [^. We assume that the both coordinates are compactified 
onto the circles with radii above. Now our claim is that the tori in the Ruijsenaars model and in branes are 
related by the T duality transformation in the both directions. Therefore, the radii in the Ruijsenaars torus 
are a R^,a g^. According to the general rules, the eigenvalues of the monodromies around the cycles are 
transformed into the positions of n branes along the corresponding directions. To be precise, now we have n 
5D branes of the TIB theory with the (worldvolume) a;®-positions which are localized at pi in the x^ direction 
and n Dl branes with x^ worldvolume coordinates which are localized at Xi along x^ . It is assumed that 
our 4+1 theory is defined on first five coordinates of the worldvolume (x", x^, x^, a;^, a;^, a;^) of n parallel D5 
branes and positions of branes pi in the x^ direction correspond to the expectation value of one adjoint real 
scalar field on its worldvolume. Let us also note that the arising configuration is the simplest example of the 
"polymer picture" from in the case when there are no "background" branes. 

Let us now switch on the interaction in the many-body problem. This means that we add the massive 
adjoint hypermultiplet in the field theory approach, add the marked point with a nontrivial monodromy around 
in the CS picture, or add "background" branes in the string language. In the CS theory, the monodromy 
relation modifies to 

where the monodromy gc. corresponds to the CP(n~l) coadjoint orbit at the marked point. After the T duality 
transformation, one gets the brane picture with a nontrivial distribution of the flux of the magnetic fleld on 
the torus surface. Actually, the flux distribution is dictated by the structure of the monodromy gc, and one 
has equal fluxes M on all plaquettes except the ones on the "diagonal" 5ij. This implies a constant magnetic 
fleld through the torus and n "vortexes" with the line structure on the diagonal which locally compensate the 
external magnetic fleld. Due to the external field the "brane polymer" is deformed and the local deformation 
on the i-th site 5pi{xij^i — Xi) is of order M in agreement with the Lax matrix structure. Note that the set 
of fluxes through the torus can be interpreted as the existence of background branes which are localized at a 
single point on a T dual torus. 



2.3 Type IIA/M theory perspective 

One of the possible deals with the Seiberg-Witten theory is to look for its origin in the M theory branes Q. 
In order to get a four-dimensional gauge theory with eight real supercharges, i.e. N = 2 SUSY, we consider 
the field theory on the world volume of the M5-brane. Its worldvolume splits into the flat J?* which is the 
space-time of the four-dimensional theory and a compact surface S of flxed genus g holomorphically embedded 
into the four-dimensional part Q4 of the transverse space. 

After compactification onto the M theory circle this single fivebrane looks as a set of Dirichlet fourbranes 
stretched between NS5-branes. 

Let us briefiy remind the basic ideas of the geometrical approach to the Seiberg-Witten theories ]17|, 
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po| . In order to get N — 2 four-dimensional SU{n) gauge theory, one has to consider the Type IIA string 
compactification on the Calabi-Yau threefold near the An-i singularity. Blowing up the singularity results in 
changing the complex structure of the mirror IIB model: 

xy^P^iw) (2.5) 

If we have more than one singularity - that is, more than one gauge group, the appropriate picture involves 
an ^p-singularity on the base. Gradually blowing it up, we reduce its rank by one: 

+ P„^(X)w^-'' + xy = (2.6) 

At the very end, one gets the "product of gauge groups" : 

+ (\)wP-' + P„, (X)w^-^ ... + xy = Q (2.7) 

Incorporation of the fundamental matter means taking coupling of the SUijif) group to zero, i.e. the size of 
the two-cycle to infinity. To have fundamental matter for each of the SU{n) factors, replace the base P^ by 
the chain of trivalent vertices, with fibers being the chains A-m-i of "gauge" P^'s and the chain Am^-i of 
matter P^'s over each "gauge" sphere. 

It coincides with the Type IIA picture of Witten's construction described above. To ensure they are 
the same theories, let us briefly revise the appearance of a fivebrane in geometrical approach. According to 
p^ , the ^n-i-singularity on the fiber over P^ in the Type IIB string theory is T-dual to the n coinciding 
fivebranes in the Type IIA theory. We can regard our ADE-type singularity as a singularity in the local model 
of K3. Then ( ^.5| ) in each root of P„{w) describes a degenerated torus. The monodromy r ^ t + 1 around 
such singular fiber corresponds to p ^ p -|- 1, i.e. B ^ B + 2-k from the dual IIA point of view [^, In 
order to see that this region carries the unit H charge, integrate it over a three-cycle consisting of a circle 
around the point with degenerated fiber times the fiber T^: 

— / H = — = lsi= 1 (2.8) 

2^ Jt^xSI Jt^xSI 27r 

It has a straightforward generalization useful for the group product case that each singular point on the 
base corresponds to a Ai'S'5-brane in the brane language These singularities (fivebranes) have also clear 
interpretation in integrable models. Their counterparts are just marked points on the bare spectral curve. 
Although the elliptic models are not the subject of the present paper, they give an insight of this connection. 
Consider the class of Hitchin systems defined on torus as the bare spectral curve: (spin) Calogero model, 
(spin) Gaudin system, etc. Our conjecture is that all these models have a corresponding brane interpretation 
with marked points corresponding to the Ai'S'5-brane. 

Another picture is actually formulated by Witten (see also ^^) who suggested to treat all branes in 
the IIA configuration as projections of a single M5 brane wrapped around the full spectral curve (not bare 
one as in the IIB picture) via different sections. 

Let us remark the only notion of the spectral curve is not sufficient to talk of integrability. Actually, 
we need an integrable dynamics that can be realized on a set of points (n for GL{n) theory) moving on the 
spectral curve. In fact, it is right their dynamics that is integrable and is linearized on the Jacobian and just 
these degrees of freedom have been missed so far. 

Thus, in order to get a complete integrable system we need the new ingredient - a set of points on the 
spectral curve which can be represented by the additional KK excitations. In lOd IIA picture they correspond 
to the nonperturbative corrections which are interpreted as DO branes on D4 branes. But of course it is not 
a whole story - one has to introduce Kqcd via a set of "magnetic fluxes". The integrable system we are 
trying to get in the brane terms and which governs the answer in the N = 2 theory without matter is the 
periodic Toda chain. It is convenient to start with the theory with adjoint matter and send the mass of matter 
to infinity. Then, due to the dimensional transmutation phenomena a new scale - Aqcd = A/ exp(— ppgj) 
emerges. 
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In the IIA brane picture, we have n D4 branes stretched between two NS5 branes. These D4 branes 
provide the worldvolume for d=4 theory. The a;* + ia;^-coordinates of D4 perturbatively provide the momenta 
in the integrable system. There are also n DO branes on D4 branes which would be ordered along the 
direction. In the Toda case, their coordinates can be naturally identified with 

= jA + 0, (2.9) 

where A is a constant which defines the equilibrium position, (f>i represent the small fluctuations corresponding 
to the Toda coordinates and the variables Xi can be considered as the coordinates in the Calogero system 
(which has no any "lattice structure", that is to say, ordering). Note also that, by comparing to the field 
theory data, one gets nA — g-i^j^j^ making our interpretation along x'' direction natural. DO on D4 behaves 
like a monopole, so we obtain something like the "monopole ring" which lies on the "diagonal". 

The possible explanation of the local Toda Lax operator now looks as follows. The element Ln tells that 
the DO is placed on the D4 brane so we can assign the same a;* + ix^ coordinate to it. Actually the whole 
issue of the Lax operators concerns the structure of the fermionic zero modes in the topologically nontrivial 
background, so the very possibility of 2 x 2 representation tells that we have no more then 2 fermionic zero 
modes here. The elements 1/12,^21 correspond to the nontrivial matrix elements of the spectral parameter 
(2.3) between the fermionic zero modes localized on the neighbor DO branes. At last, we would conjecture that 
vanishing of L22 reflects the fact that only one chiral zero mode exists for a monopole under consideration. 

In what follows we exploit this monopole-related interpretation of the brane picture. Let us emphasize 
that we would fairly talk about the monopoles in the SU{n) theory not about n monopoles of the SU{2) one, 
the latter seems to appear because of some duality. Later we will confirm that just these SU{n) monopole 
system can be immediately generalized in a more complicated situation. Let us also remark that evidently we 
have no the 2x2 monopole picture for the Calogero system which is in a perfect agreement with the absence 
of the second Lax representation (see section 5) . 

To complete the section outline once again that in any picture three main ingredients have to be included. 
First, one guesses the set of branes which have to provide the ultimate worldvolume for the d=4 theory with 
SU{n) gauge group - these are n D7 in the IIB/F or M5 in the IIA/M pictures. Second, one has to add 
"instantons" which are M2 and KK excitations in the M theory and finally the regulator degrees of freedom 
are introduced - these are Wilson line on the torus in IIB/F or some "magnetic flux" analogues in M theory. 



3 Simple gauge group and SL{2) spin chains 

In the section we discuss integrable systems that correspond to some N = 2 SUSY YM theories. Among these 
theories there are a few which are UV-finite. Integrable systems corresponding to the UV finite theories softly 
perturbed by mass terms serve as "reference systems" for us, since all other (asymptotically free) theories can 
be obtained from the UV-finite ones by some degenerations. 

The defining property of integrable model behind UV-finite theory is that it is associated with elliptic 
"bare curve" whose modulus r plays the role of the bare coupling constant in the theory. In some cases 
(theory with adjoint matter hypermultiplet Q) this bare curve is an underlying Riemann surface where the 
Lax operator is defined. In other cases, considered in this paper, some coefficients (coupling constants) turn 
to be some modular forms of r. In this latter case, one can ignore the bare curve in the four-dimensional 
consideration, just keeping in mind that it lives in the 11-12 dimensions of F theory and dealing with the 
coupling constants as just numbers (that are properly adjusted when doing the dimensional transmutation). 

We start with the simplest case of the theory with the SU{n) gauge group. The integrable structures 
corresponding to these case are well-known [^, ^, ^ and we reproduce it here merely in order to have some 
basic simple example. 

This theory has zero /3-function if one adds 2n massive matter hypermultiplets in the first fundamental 
representation. Before going into details let us note that in the brane picture B this case corresponds to 
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the two parallel NS5-branes with n D4-branes stretched between them and 2n additional D6-branes placed 
between NS5 that provides the UV-finiteness of the theory. In accordance with our general rules, it corresponds 
to the "reference" integrable system. This integrable system is the periodic inhomogeneous SL{2) XXX spin 
chain with period n has to be slightly twisted to make possible the dimensional transmutation to the 
asymptotically free theories. Let us note that, although throughout the paper we use the term "inhomogeneous 
SL(p) magnet", it is equivalent, for the rational (XXX) case, to the homogeneous GL{p) magnet. 

3.1 UV-finite theory 

The periodic inhomogeneous sl{2) XXX chain of length n is given by the 2x2 Lax matrices 

3 

U{X) = {\ + X,)■l + J2Sa,^■CT'' (3.1) 

a = l 

a"' being the standard Pauli matrices and Ai being the chain inhomogeneities, and periodic boundary condi- 
tions. The linear problem in the spin chain has the following form 

L,(A)*.(A) = *.+i(A) (3.2) 

where ^'i(A) is the two-component Baker- Akhiezer function. The periodic boundary conditions are easily 
formulated in terms of the Baker- Akhiezer function and read as 

*»+n(A) = ^w'i'.iX) (3.3) 

where w is a free parameter (diagonal matrix). One can introduce the transfer matrix shifting i to i + n 

T(A) = L„(A)...Li(A) (3.4) 

which provides the spectral curve equation 

det(T(A) + w ■ 1) = (3.5) 
and generates a complete set of integrals of motion. 

Integrability of the spin chain follows from quadratic r-matrix relations (see, e.g. |]l^ ) 

{L,(A) ®L,(A')} =<5,, [r-(A-A'), L,(A) ® L,(A')] (3.6) 

with the rational r-matrix 

-I 3 

a = l 



The crucial property of this relation is that it is multiplicative and any product like (3.4) satisfies the 
same relation 

{r(A) ? T(A')} = [r(A - A'), r(A) ® r(A')] (3.8) 
The Poisson brackets of the dynamical variables Sa, a — 1,2,3 (taking values in the algebra of functions) 



are implied by ( 3^ ) and are just 

{Sa, Sb} ^ -i^abcSc (3.9) 

i.e. {Sa} plays the role of angular momentum ("classical spin") giving the name "spin-chains" to the whole 



class of systems. Algebra (3.9) has an obvious Casimir function (an invariant, which Poisson commutes with 
all the spins Sa), 



3 



74-2 =8" = ^5'aS, (3.10) 



The spectral curve (3^) can be presented in more explicit terms: 

^2 ^ ^^j^^x)w + detr(A) = (3.11) 
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In accordance with ([4. j), the last term defines masses of the hypermultiplets. Since 



one gets 



det2x2L»(A) = [X + X.f-K'' 

n 71 

det2x2T(A) = J|det2x2Li(A) = H ((^ + ~ ^0 = 

n 

= n(A-m+)(A-mr) 

where we assumed that the values of spin K can be different at difi^erent sites of the chain, and 

mf = -Ai ± Ki. 



(3.12) 



(3.13) 



(3.14) 



While the determinant of monodromy matrix ( 3.12 ) depends on dynamical variables only through Casimirs 
Ki of the Poisson algebra, the dependence of the trace Tr2x2r(A) is less trivial. Still, it depends on S'i*' only 
through Hamiltonians of the spin chain (which are not Casimirs but Poisson-commute with each other) - see 
further details in 

Let us note that we have some additional freedom in the definition of the spin chain and the spectral 
curve. Namely, note that r-matrix (3.7) is proportional to the permutation operator P[X ®Y) = Y ® X . 
Therefore, it commutes with any matrix of the form U ® U . Thus, one can multiply Lax operator of the 
spin chain by arbitrary constant matrix without changing the commutation relations and conservation laws. 
Moreover, we can also insert a constant (external magnetic field) matrix V into the end of the chain (into the 
n-th site). This is the same as to consider more general boundary conditions - those with arbitrary matrix 

. This is why such a model is called twisted. 

The described freedom allows one to fit easily the form of the spectral curve proposed in ^ 

.-^^P(A), 



2h 

TJli > ] 

n ^ 



h{T)^ 



where r is the bare curve modular parameter and Qi are the theta-constants [ ^ . 
It can be done, e.g., by choosing matrices U and V to be^ 




I 



V ■ 



n- 



h{h+i) 





(3.15) 



(3.16) 



/ 



n- 



h{h+i) 
-1 



(3.17) 



Let us remark that further extension of the theories with the fundamental matter would lead to the 
(twisted) XXZ and XYZ chains for the 5d and 6d theories correspondingly. 



3.2 Pure gauge theory 

Now let us briefly consider the degeneration of our reference system to the Uf < 2n case. This can be done in 
the standard way Q| by taking I masses mi, mi to infinity while keeping Aqqjj = e"^^mi...m; finite. After 
this procedure the modular forms disappear so that Aqcd emerges instead. 



'To fit (3.15), we also need to shift Ai Ai — ^ mi 
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Degenerations of the system can be studied at a single site (for the sake of brevity, we omit the index 
of the site). Let us consider the Lax operator L = UL (see (3.1), (3.16)) with spins satisfying the Poisson 



brackets (3.9). We are going to send a to zero. Along this line of reasoning, there are two possibilities to get 




nontrivial Lax operator. The first possibility, when the both masses (3.14) disappear and one reaches the pure 
gauge theory, is described by the periodic Toda chain In order to get it, one needs to redefine S+ -^3+, 
then take a to zero and remove, after this, the inhomogeneity by the shift of So. This brings us to the Lax 
operator of the form (we introduce new notations So ~ S3, S± — Si i: iS2) 



(3.18) 

(3.19) 
(3.20) 

(3.21) 

(3.22) 



so that the Poisson brackets are 

{S±,So} = ±S±, {S+,S-} = 
This algebra is realized in new (Heisenberg) variables p and q 

S± = ±e±^ So=p, {p,q} = l 
This leads us finally to the Toda chain Lax operator |l2| 

^ A + p e" \ 



LToda 







and the spectral curve (3.11) takes the form 



w + TrT{\)w +1 = 



Now let us return to the second possibility of the asymmetric degeneration, when one of the masses (3.14) 



remains in the spectrum while the second one goes to infinity. One can understand from ( 3.14 ) that, in contrast 
to the Toda case, this degeneration requires a special fine tuning of the Casimir function and inhomogeneity, 
so that both of them go to infinity but their sum (difference) remains finite. In the Lax operator terms it 
can be done in the following way. Let us redefine S+ ^5+ and So — > ^So. This means that the Poisson 
brackets take the form 

{S+,S-}^-2So, {5'±,So} = (3.23) 



Now in order to preserve finite Lax operator (3.1), one needs to take care of its element Lii{X). This can be 
done by the rescaling Xi —Xi and fixing Xi + So to be a • so. This brings us to the Lax operator 



L(A) 



A + so S_ \ / A + so S- \ 
\ S+ Ai — So / \ S+ — 2So / 
The determinant of this Lax operator is equal to (A — m) where m is the finite mass 

m = 2soSo + S+S- 



(3.24) 



(3.25) 



in perfect agreement with (3.14), Let us note that 2soSo + S+S_ is also the Casimir function of the algebra 



(|3.23| ), since {S±,so} = ±S±. 

In fact, the Lax operator ( [3.24| ) can be rewritten in the form 

A + s+s_ — m s_ 
s+ 1 



L(A) = 



(3.26) 



where we rescaled the spins and used the explicit form of fixed the Casimir function 2soSo + S+S-. This 
Lax operator defines so-called novel hierarchy and is gauge equivalent to the discrete AKNS ^] and 
relativistic Toda chain |^^. This latter correspondence looks especially unexpected, since the same relativistic 
Toda chain describes the 5d pure gauge SUSY theories [p3[. 
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Let us note that another (equivalent) way to count all possible degenerations p3 is to consider 2x2 
Lax operator of the most general form linear in A, which satisfies the Poisson brackets (p.q) with the rational 



r-matrix (3.7) and determine all Casimir functions with respect to this Poisson brackets. Then, all possible 
degenerations are in one to one correspondence with the vanishing of Casimirs^. 

To complete this section, let us say a few words of the brane interpretation of the Toda variables pi. First, 
restore the dependence on Aqcd in the Toda Lax operator. It can be done in many different ways, we choose 
it to be 

/ A+p AlcDe" \ 
LToda = (3.27) 



thus the spectral curve is of the form 

+ TrT(A) + A|"ci5 = (3.28) 

Now taking the perturbative limit, i.e. putting Aqcd 0, we get that the spectral curve degenerates to the 
two sphere touching each other in the n points that are zeroes of the polynomial Trr(A). These zeroes are 
equal to pi in the perturbative limit and describe the positions of the D4-brane ends in the (x* + ia;^)-plane, 
in accordance with the picture of M. 



There is another class of the spin chain variables, namely, those describing masses of the sixbranes (3.14) 



which can also allows very rough but transparent brane interpretation. Indeed, in accordance with (3.14), one 
can associate with each D4-brane (i.e. site of the spin chain) a pair of the sixbranes, the distance between 
them being measured by the value of the Casimir function K. At the same time, the (x* + ja;^)-coordinate 
of the center of masses of these two sixbranes is naturally measured by the inhomogeneity A^. This looks 
especially reasonable, since just the inhomogeneity corresponds to the decoupling U{1) factor of the GL{2) 
magnet.^ 

Note that, within this interpretation, the procedure of the degeneration gets its pictorial explanation. It 
just corresponds to removing one of the sixbranes to infinity. To keep the other branes at finite distances, one 
needs to tune parameters \i and K so that one of the two quantities ~\i ± K remains finite. 



4 The group product case 
4.1 Reference system 

Consider now more complicated case of the gauge group which is the product of some SU{n) factors: G = 
X^i-\ SUijii). One can also add p — 2 massive hypermultiplets in the representations (ni,ni+i) and also by 
some hypermultiplets in the first fundamental representations. This case was first considered by E.Witten 
in Q and possesses a new important feature. To start with, let us note that in this case for many different 
choices of rn the suitable number of matter hypermultiplets can be added to make /3-function zero. Actually, 
one needs to add 

= 2ni - rii-i — rii+i, no = rip = (4.1) 



^We are grateful to S.Kharchev for the useful discussions on this point. 



''One more interpretation of ( 3.14 ) might look as follows. We can start with a free theory and would 
expect the monopoles in "the momentum space" by the following reason. For the free particle, there is 
the standard dispersion relation E — ±(p*^ + m^)^ which tells us that there is the level crossing on the 
multi-dimensional surface p = ±im in the complex momenta space. Would be this surface one-dimensional, 
according to the general philosophy of the Berry phase, this means that one needs to introduce two monopoles 



in "the momentum space" with their centers located at points p = ±jm. This resembles formula (3.14) without 
inhomogeneities. Furthermore, one can switch on the interaction with the Higgs field that results into the shift 
of the fermion mass due to the Yukawa coupling. Since inhomogeneities are naturally related to the Higgs 



vev's, one finally arrives at formula (3.14) 
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Figure 2: Brane picture, corresponding to the 5/(3) magnet 

hypermultiplets [^. This always can be done unless some m — m+i are too large. However, among all 
the integrable systems describing these numerous UV-finite theories there is a unique one that plays the 
role of the reference system, since all others can be obtained by degenerations. This theory is specified 
by maximal possible number of the matter hypermultiplets at given p and ni and characterized by the set 
{rii} — n, 2n, 3n, . . . ,{p — l)n. The corresponding number of the hypermultiplets is ni + Up^i — np. 

Now come to the brane picture. In this case we have p parallel NS5-branes with ni D4-branes stretched 
between the first pair of them, 712 D4-branes stretched between the second pair of the NS5-branes etc, with 
the corresponding numbers di of D6-branes. The reference theory in this language corresponds to Fig. 2. 

Its degeneration to theories with fundamental matter (D6-branes) in brane languages corresponds that 
A = + ix^ positions of several fourbranes coincide. Namely, if k of rip-i D4-branes have the same A 
coordinates as fc — 1 of np_2 and so on till the only D4-brane of Up-^k set, we can attach a fictitious D6-brane 
to the right infinite end of the D4-set with this value of A = a;* + ix^, and move it through fc — 1 fivebranes 
to the left consequently annihilating the whole set of such D4's. This mechanism will be explained in details 
in section 5 via brane creation phenomenon. 

In further assigning brane configurations with integrable systems, it will be of importance that, with any 
described branes system, there can be associated the spectral curve of the form |6| 



p- 



(4.2) 



gn2—di (A)Ji(A)w''-' + ...4 

k 

Here (/i(A) is a polynomial of degree i depending on the Coulomb moduli^ 
and hypermultiplet masses, and 

n (^-"^i) (4-3) 

are pure mass terms (mj 's are masses of the hypermultiplets in the fundamental representations) . 

Let us note that this spectral curve gives another justification for our choice of reference system. Indeed, 
let us consider "the pure gauge case' 'Qwhen aU Ji(A) = 1 g. Then, only the moduh space of the reference 
spectral curves contains the Z„-symmetric curve, and, in this sense, is maximally possible complete (although 



^Hereafter, we call "Coulomb moduli" the moduli corresponding to the vev's of the fields in the vector 
(gauge) multiplet. 

^To simplify further terminology, hereafter we call the system that contains no fundamental hypermultiplets 
"the pure gauge theory" , despite there are still presented p—2 hypermultiplets in the representations (n^ , rii+i) . 
Our terminology refers to the fact of absence of any D6-branes and semi-infinite branes and takes its origin 
in the SU{n) case. 



12 



the moduli space of curves that correspond to an integrable system is always not complete for sufficiently large 
genera and covers typically dimensional subspace of the complete {3g — 3)-dimensional moduli space; the g 
moduli are just the Hamiltonians of integrable system and, therefore, the dimension of the configuration space 
is exactly g). 

To clarify this point, let us introduce the new variable Y — w — Then, only with the choice nt = kn 
one can match the polynomials gt so that the curve takes the form 

Y" = R{X) (4.4) 

where R{X) is a polynomial of degree np, i.e. this curve is of genus 

ff=f^-l)(P-l) (4-5) 



2 

The same genus can be obtained immediately from the brane picture, if one implies that the D4-branes are 
blown up to the handles of the Riemann surface, while NS5-branes are the spheres that are glued by these 
handles. Then, n D4-branes between first two NS5-branes give n — 1 wholes, 2n next branes give 2n — 1 holes 
etc. - totally, it is exactly — l) (p — 1) holes. 

In fact, this interpretation of the Riemann surface (spectral curve) as the brane graph implied in is 
justified by the calculation of genus of the general curve 

+ g„, {\)w^-^ + (A)^^*'-^ + ... + ffn, {\)w^-'' + ... + (A) = (4.6) 

The degrees of polynomials in this expression are restricted by the conditions m < n2 < ... < n-p-i and 
ni > n2 — ni > na — n2 > ... > rip — np_i. The second condition is the requirement for all the /3-functions 
to be non-positive, i.e. for the theory to be asymptotically free or UV-finite. These conditions allow one to 
calculate easily the genus of the curve. For doing this, we use the elegant formula (see, for example, [Q)|^ for 
the genus of the non-singular projective curve given in the affine map by the polynomial F{w, A) = ^ aijW^X' . 
The genus is merely equal to the number of the integer points in the Newton polygon A for F, i.e. the convex 
shell of the points (i, j) G Z such that aij 7^ - see Fig.S.p] 

Since powers of w represent resolution of singularity on the base, and the degree of each polynomial g^. 
corresponds to resolution of A^. _i singularity on the fiber, by construction, A is nothing but toric polyhedron 
responsible for the local geometry of Calabi-Yau threefold [Q. 

Following each interior node of A define a compact divisor, thus generating a homology group. As 
it was explained in details in |2^, Riemann surface E faithfully represents all the data of Calabi-Yau cycles, 
with correspondence between cycles and forms. In other words this means that each compact divisor counts 
an element of i/i(E). 

In the present paper, from the integrable system point of view, and in from the geometrical side, it 
was independently inferred that to have UV-free four-dimensional theory, i.e. each di > 2ni — rn-i — rn+i, 
means to consider only convex A. 

Using the conditions onto the set of {rn} (the second condition just guarantees the convexity), one can 



easily get that the genus of the curve (4.6) is equal to 



p-i 

^ n. - p -M (4.8) 



'^"We are grateful to A.Chervov who has drawn our attention to this formula. 

'^^This formula is proved by noting that, for each pair inside the polygon, there is the holomorphic 

differential 
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Figure 3: Newton polygon A 



Let us note that this result coincides with one calculated from the brane picture as described above, since 
m D4-branes between first pair of the NS5-branes gives m — 1 holes etc. Remark also that in accordance with 
expectations, the degree of the last term does not enter the final answer, which means that the genus is not 
effected by the D6-branes (matter hypermultiplets) added. 

Now, using the formula ([f .8|), we can determine the reference system as corresponding to the curve of 
the maximal possible genus at given p and ni — This genus is equal, in accordance with (4.8), to 

(p-1). 

Below we describe the reference integrable system that is the periodic inhomogeneous SL{p) (or homo- 
geneous GL{p)) spin chain with period n. It is described by the curve (IS) with all but Jp-i equal to 
1. 



4.2 UV-finite theory 

The inhomogeneous SL(p) spin chain is similar to the SL(2) spin chain considered above and given by the 
p X p Lax operator 

L^{X)=K''''Sa.iXa+{\ + X^)■l (4.9) 

at the i-th site. Here Sa,i are dynamical variables^ Xa axe generators of the SL{p) algebra and K°'^ is its 
Killing form. We always choose below Xa to lie in the first fundamental representation of SL{p). In complete 
analogy with (3^), this Lax operator shifts the Baker- Akhiezer function to the neighbor site: 

L.(A)*.(A) = *,+i(A) (4.10) 

and periodic boundary conditions read as 

*»+„(A) = -w<f,{X) (4.11) 
Introducing the transfer matrix shifting i to i + n 

T(A) = L„(A)...Li(A) (4.12) 
and using boundary conditions, one gets the spectral curve equation 

det(T(A) +u> ■ 1) = (4.13) 



'^^Note that, by the shift of w, (?ni(A) can be cancelled. For the reference curve, however, this shift does not 
change the degrees of all the rest polynomials gni(A). 

^^To reproduce the SL{2) case considered above, one needs to make the replacement So 2So- 
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In fact, (4.13) is a polynomial in w and A with the coefficients that give an ample set of integrals of motion. 



The Poisson structure is in charge of integrability of the system and is given by (3.6) with the rational 
r-matrix 

,(A) = ^^!!^^i®^ = ^ (4.14) 
A A 

where P is the exchange operator P{X igiY) = Y ® X and the last equality is correct, since we consider the 
fundamental representation. This gives rise to the Poisson brackets of Si 

{Sa,Sb} = -f:tSc (4.15) 

where /^j are the structure constants of the SL{p) algebra. The Poisson structure is, generally speaking, 
degenerated so that its annulator is the algebra of Casimir functions that are defined to be invariant with 
respect to the coadjoint action. Therefore, one needs to restrict the Poisson bracket to the orbits of this action. 

Traces of degrees of Lax operators are Hamiltonians which mutually commute with respect to this Poisson 
structure. It can be easily understood since ( |3.6| ) is multiplicative and, therefore, products of Lax operators 
satisfy the same Poisson brackets. 

Now one can use the explicit formulas for the Lax operator to obtain the following spectral curve 

w'' + gi{X)w''-^ + ... + gk{\)w^-'' + ... + 3p(A) = (4.16) 

where gk is a polynomial of degree kn. Manifestly, gi{\) = Trr(A), gk ~ detMj''' (A), gp{\) = 
detT(A), where Mj''\\) is the matrix obtained from the transfer matrix by removing k columns and k rows 



that intersect at the diagonal elements given by the multi-index I. This spectral curve is exactly (O) with 
all but Jp-i put equal to 1. 

Apart from the spectral curve, there is another key ingredient of the integrable system which is of great 
importance in the supersymmetric theories. This is the meromorphic differential dS ^ |35j that generates 
the spectrum of the BPS states in the theory [Q, |^ . The defining property of dS is that its variations over 
Coulomb moduli gives holomorphic 1-forms. In the spin chains, this differential has the form 

dS = A— (4.17) 

w 

Let us check that this is, indeed, the proper differential, i.e. its variations over the Coulomb moduli give 



the holomorphic differentials. Consider the curve (4.16) (or ([l.q)) and write it in the form F{w,X) = 



'^ij^^^'' ~ 0- The differential dS depends on moduli aij of the curve only through the dependence of 
A determined by the equation F{w, A) = 0. Thus, we get 

ddS _ d\ dw _ OF dw _ w'-^X^dw 
duij duij w daij wF'^ F( 

Now it remains to note that the coefficients ap-k.n^ ^■■'^ ju'^t fixed numbers (unities)|^ that are not to be 
regarded as moduli, while gp{\) — J^(A — rrii) does not contain the Coulomb moduli and only masses. Then, 
using (^^) we finally get that the variations of dS over the Coulomb moduli do really generate all the 
holomorphic differentials. 

In order to check that the inhomogeneous periodic SL{p) spin chain is really the reference system we 



are looking for, one needs to check that it properly degenerates to reproduce the curve (4.2) with arbitrary 



Ji{X). Indeed, from (4.16) one reads off that Jp-i = detT(A) — ]^"detLi(A). Each determinant detLi(A) 
is a polynomial of degree p with coefficients being Casimir functions. On the generic orbit, all zeroes of this 
polynomial and, therefore, masses are different. Thus, we obtain Jp_i = ri"nr('^ ~ mi.r), where mi,r refers 
to p different roots of detLi{X). 

In order to get generic J;, one needs to consider some coinciding masses at the same site. If only two 
masses coincide, one gets only Jp-i ^ 1 and Jp-2 7^ 1. This case corresponds to a specific set of orbits with 



This property of (4.16) follows from the manifest form of the Lax operator (4.9) 
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one relation between Casimir functions. Since there are p — 1 (= rank SL{p)) independent Casimir functions, 
there can exist maximally p — 1 relation between them. The corresponding special orbits describe Ji(A). 
Certainly, at different sites, one can consider different degenerations, giving (?p(A) the general form as in (4.2). 



Now this is nontrivial test of the whole construction to check that the spectral curve (4.16) with all these 
degenerations gets the form ( [4.2| ), i.e. as soon as two roots of the polynomial gp{\) coincide, the polynomial 
gp-i{X) gets simple zero at the same point; as soon as three roots of the polynomial gp{\) coincide, the 
polynomial gp-i{X) gets double zero at the same point, while the polynomial gp-2 gets simple zero at this 
point etc. This is, indeed, the case and in the Appendix we demonstrate this by the direct calculation in the 
simplest case of the SL{3) spin chain. 

To conclude this subsection, let us note that the interpretation of the mass formulas for the sixbranes 
proposed in sect. 3 is naturally continued to the SL{p) magnet. Now a site of the spin chain should be associated 
with one D4-brane suspended between the first pair of the NS5-branes, two D4-branes suspended between the 
second pair, p — 1 D4-branes suspended between the last pair. With this last set of the p — 1 D4-branes 
is then associated p sixbranes, with their coordinates being described by the masses. In this situation, again, 
the inhomogeneity describes the (a:"* + ia:®)-coordinate of the center of masses of the sixbranes, and bringing 
one of the branes to infinity describes the procedure of the degeneration. 



4.3 Pure gauge theory 



Now let us discuss decoupling of the fundamental hypermultiplets, that is to say, degeneration of the spin 
chain. As before, the degeneration can be manifestly done using the freedom of multiplying the Lax operator 



by any constant matrix, since the r-matrix (4.14) is again the permutation operator and commutes with 
the matrix U x U. This procedure looks literally the same as in the 5*1/(2) case. Thus, let us now discuss 
only the maximally degenerated systems corresponding to the pure gauge theories. These axe described by 
5L(p)-generalizations of the Toda chain system. 

As in the SL{2) case, it is sufficient to consider the Lax operator at a single site. In the SL(p) case there 



are p — 1 different possible degenerations of (4^) generalizing Toda system and p — 1 different Lax operators. 
They contain the spectral parameter A in the only one diagonal term, in two diagonal terms etc. till p — 1 
diagonal terms and have unit determinant. Varying this different Lax operators along the chain, one can 



reproduce polynomials of different degrees in (4.2) with all A (A) = 1 (i.e. all di — 0), the degrees being 
subject to the inequalities ni < 712 < 



< Up-i and ni > n2 — ni > — n2 > ... > Up 



This 



really corresponds to the pure gauge theory, see . These different possible polynomials correspond to several 
possible partitions of D4-branes between background D6-branes. 



As we already showed in s.3.2 (very explicit formulas for the 51/(3) case can be found in Appendix), one 



can obtain the different degenerations just multiplying the Lax operator (4.9) by constant diagonal matrices 
U with 1, 2, p — 1 non-unit entries equal to a, rescaling some spin operators S — » and then taking a 
to zero. Here we demonstrate the equivalent way of doing briefly mentioned in the very end of s.3.2 [psj. 



Namely, let us start with the general p x p Lax operator that satisfies the Poisson brackets (3.6) with the 



rational r-matrix (4.14) 



{Uj{\),Lki{y)} 



A- 



{Lk](\)Ui(p) - L«(A)Lfcj(/i)) 



(4.19) 



Consider specific representations that are described by the Lax operators of the form (we consider, for a 
moment, the homogeneous GL{jp) system instead of inhomogeneous SL{p) one) 



L,,(A) = a,A ■ l-f^,,, TrA = 



These Lax operators can be produced from (4.£) multiplication by a constant matrix. From 
read off the Poisson brackets for the matrix A 



(4.20) 
1£|) one can 



{Aij , Ajjj} — akAiiSkj ~ aiAkjSii 



(4.21) 



16 



with tti being Casimir functions. From these Poisson brackets one concludes that, for an element Aij to be a 
Casimir function, there should be fulfilled the condition = a, = 0. Now let us look at the determinant of 



the Lax operator (4.20). As was explained above it has the form 

detL(A) = aia2...ap\^ + CiX^^'^ + CaA''"^ + ... + Cp (4.22) 
where d are the Casimir functions (this formula specifies their choice). In the general SL{p) magnet, all 



at 7^ and can be put equal to 1. In this case, the determinant (i.22) is a general polynomial of degree p with 
p zeroes. Degenerating of this system implies that some (say, k) of these zeroes (masses of the hypermultiplets) 
are taken to infinity, i.e. some = and some (fc) first Casimir functions vanishes too. The pure gauge 



theory corresponds to the limiting degeneration when the determinant (i.22) is equal to constant (unity). 
This is reached by putting some a; and all Casimir functions but the last equal to zero. Let us note that there 
exist exactly as many as p — 1 classes, since one can put 1, 2, p ~ 1 a^'s equal to zero. This completes the 
construction. 



5 Two Lax representations 

5.1 2x2 versus n x n representations 

So far we dealt with "lattice" Lax representations for the SL{p) groups formulated in terms of p x p 
matrices whose product along the chain - transfer-matrix - gives rise to the spectral curve. However, it is 
often needed another type of the Lax representations given by n x n matrices that is adequate for the (spin) 
Calogero (Toda) systems (and, more generally, Hitchin systems [ p^ , p^). In this case, the Lax operator 
itself generates the spectral curve, while no analog of the transfer-matrix is known. 

However, it sometimes happens that the same integrable system can possess both Lax representations 
resulting in the same spectral curve and different but equivalent phase space descriptions for the system. In 
the brane language it can be formulated as the equivalence after the rotation of two brane configurations. An 
example of this situation will be considered in the next section. It appears that this issue can be formulated 
in a relatively rigorous way. 

In this and partly the next sections, we discuss interrelations between different Lax representations and 
their interpretation in brane terms. 

The simplest relevant example is the Toda chain which has the 2x2 Lax representation (see s.Toda) 



LToda 



\ + P 







(5.1) 



with the linear problem (3.2) and the boundary conditions (3.3). This system can be reformulated in terms of 



n X n matrix. Indeed, consider the two-component Baker- Akhiezer function 



problem (3.2) can be rewritten as 



Then the linear 



(5.2) 



and, along with the periodic boundary conditions (B.3) reduces to the linear problem £(ui)"I> — A"I> for the 
n X n Lax operator 

w 
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Figure 4: Different representations of fundamental matter first by inserting sixbranes, and 
then transformation of them into semi-infinite fourbranes 



with the n-component Baker- Akhiezer function $ = {e '"'^^i/'„}. This leads us to the spectral curve 



det„ 



, {C{w) - A) = 



(5.4) 



which is still equivalent to the spectral curve (3.5) 



The symmetry of two different representations of Lax operator by 2 x 2 or A'^ x A'^ matrices can be clearly 
reinterpreted in string theory language. 

From the brane pictures point of view the two Lax representations are nothing but the choice of parametri- 
zation of the fivebrane worldvolume E. Of two holomorphic coordinates (which are Cartan elements of scalars) 
A = + ix^ and s = + ia;^", |^ we can regard any we wish as a spectral parameter. In the "standard" 2x2 
representation, w = exp ^ is a spectral parameter, and the spin chain is the appropriate integrable system. 
The corresponding typical brane picture is depicted in Fig. 4. Another possibility is to choose A as a spectral 
parameter. It leads to the n x n representation. 

Geometrical image of this phenomenon is the base-fiber symmetry of the Calabi-Yau threefold . As it 
was already mentioned earlier, in the "standard" 2x2 language w parametrizes the base Am singularity and 
counts the number m of gauge factors. A is a coordinate in the fiber. If each coefficient at w* is a polynomial 
of degree n, we deal with the SU{n)"^ gauge theory with n extra fundamentals at each end of the chain of the 
gauge groups. The duality SU{n)™' ^ SU{m + l)^~^ proposed in 0] is nothing but the m-|-lxm-|-l nxn 
symmetry of the Lax representations. Under this symmetry, the Newton polygon (Fig. 3) refiects so that axis 
i interchanges with axis j. 



Note that the holomorphic differential 



d5 = A^ 

w 



\ds '• 



-sd\ 



(5.5) 



also suggests this symmetry. 



5.2 Sixbranes versus semi-infinite branes 

Below we argue that the representation of the fundamental matter by the D6 branes and the semi-infinite D4 
brane has difi'erent interpretation in terms of integrable systems. Namely, we would like to propose that the 
D6-branes are natural in the 2x2 representation, while the semi-infinite D4 branes in the nxn one. 

We start with making use of the effect (first found by Hanany and Witten [^) of creating a new brane 
when two other branes pass through each other. 

In the paper, the authors considered the three dimensional theory in the common longitudinal directions 
of the branes: 

^^Equivalently, s = a;^ -f- ix^ in the Type IIB language. 
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(5.6) 



By T-duality in the third direction, this setup gets mapped into Type IIA brane picture |6| described above 
(2.4) so that the statement of creating a new D3-brane when D5-brane passes through NS5-brane transforms 
into the following one: moving D6-brane through NS5-brane creates a new D4-brane stretched between them. 

As it was also mentioned in js^, the a;''-coordinate of the sixbranes apparently does not play any role, 
at least, in the low-energy physics we are interested in. Therefore, one can take these D6-branes to infinity, 
stretching a new D4-brane each time we pass through a NS5-brane. Moreover, we can move every particular 
sixbrane independently to the left or to the right up to our choice. The final configuration, shown in Fig. 2, 
includes the original setup of fivebranes and fourbranes and very long fourbranes (with a sixbrane attached 
to the end of each D4-brane) instead of the sixbranes. Then, if the sixbranes (endpoints of the D4-branes) 
are quite far from the picture we are dealing with, we can replace this setup of branes by the same without 
any sixbranes, since they do not affect the physics. For example, the beta-functions of the gauge theories 
remain the same, since the presence of the D6-branes at infinity does not change the relative bending of the 
fivebranes. In other words, it means that we are going to the infinity of Q4, which is Asymptotically Locally 
Flat. 

The use of this phenomenon can be found in other columns of Fig.l. For instance, the existence of two 
different brane pictures which reveal the same low-energy physics suggests the existence of two different inte- 
grable counterparts. Really, the appropriate counterparts for the models with sixbranes are spin chains, while 
the semi-infinite fourbranes find their description in deformations of the Toda system by nontrivial bound- 
ary conditions suggested by I.Krichever and D.Phong |^^. Though at the low-energy level both systems are 
equivalent, the difference must appear in higher terms in the field theory, or higher times in the corresponding 
integrable systems. This also responds the question about the role of the a;^-coordinates of the sixbranes. 

Let us more comment the same point from a slightly different view. 



In the 2x2 representation, the spectral curve is (3.22). This spectral curve describes the pure gauge 
system obtained by degeneration of the (D6-branes -I- pure gauge) theory (spin magnet). Therefore, it is 
natural to identify the variable w in (B.22) with the variables y or z in Q. These latter can be considered as 
describing one of the complex structures of the multi-Taub-NUT solutions 



yz = 



(5.7) 



where d are (a;''-|-i2;^)-coordinates of the sixbranes. It is of crucial importance that y and z never go to infinity 
at finite A and, thus, describe the brane configuration that have only points at finite distances whenever A is 
finite. 



In contrast to the 2x2 representation, in the n x n representation the spectral curve has the form (5.4) 



w+- = P(A) 

w 



(5.8) 



and is obtained from the Lax operator (5.3) given on a cylinder (double-punctured sphere). Let us note that 
the curves (3.22) and (5.8) are identical only if w 7^ 0. If, however, one considers the projectivization of these 



spectral curves, they will become identical, with ( 3.22 ) and (5^) corresponding to different maps. From the 
point of view of the Toda theory, this merely corresponds to different choices of gradation in Lax operator. 

It is clear that the variable w lives on the cylinder with radius R which is identified with the radius of the 
compactified eleventh dimension in M theory M 



t 



(5.9) 
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The zero values of w = t describe infinitely far points of the brane configuration. If now one adds the matter 
hypermultiplets in the fundamental representation, the spectral curve takes the form 

w+^ = P{\) (5.10) 



equivalent to (3.11) unless w = 0. However, there will be the points on the curve with zero values of w — t (i.e. 
infinite values of physical coordinates) located at finite A's. These points describe the semi-infinite D4-branes 
1^ and can be thought of as infinitely long handles, or punctures on the Riemann surface that correspond to 
the matter fundamental hypermultiplets in the standard picture pm. 



The curve (5.1C) can be immediately obtained from the n x n representation by the simple deformation 



of the Lax operator (5.3) just filling up the first and the last rows of the matrix (5.3) with any constant 
(independent on A and proportional to and to w respectively) entries |^^. This Lax representation can 
be just obtained as already mentioned non-local boundary condition in the Toda chain jssj . Certainly, one can 
equally add non-zero entries to the first and the last columns. These two possibilities correspond to attaching 
the semi-infinite branes to the right or to the left NS5-branes respectively. We return to this construction 
again in the next subsection. 

In a word, we would associate 2x2 spin chains with the system with matter fundamentals realized via 
inserting the sixbranes. On the contrary, the n x n representation is associated with the realization of the 
matter hypermultiplets via the semi-infinite branes. 

There is another argument in favor of this identification that is related to the integrable system ^ brane 
configuration correspondence. Indeed, it is natural to consider the infiuence of the semi-infinite branes as a 
boundary efi^ect like the emergence of Q{X) in the paper is due to the non-local boundary conditions. At 
the same time, the sixbranes are more naturally associated with the "local" perturbations of the D4-branes 
suspended between NS5-branes, i.e. with additional degrees of freedom in the spin chain. This nicely fits the 
recent observation that the semi-infinite D4-branes can be more naturally attributed to the meson not 
quark degrees of freedom. 



5.3 n X n representation in the higher rank magnets 



In this subsection we discuss a way to obtain n x n representations by degenerating the SL{p) magnet. On 
this way, we reproduce the Lax operator of that corresponds to bringing to infinity all but the very left 
two NS-branes. Then, the resulting system looks as the SU{n) SYM theory with massive hypermultiplets 
realized via semi-infinite branes attached to the right NS-brane. 



To begin with, let us note that, starting from eq.(5.4), one can maximally reproduce the curve (5.8) with 
polynomials P{X) of degree n and (5(A) of degree 2n — 2. Thus, within the n x n representation, one can not 
reproduce the reference system (that requires the degree 2n of the polynomial Q(A)) but can approach to it 
as close as 

possible^ 

Now let us consider the pure gauge limit of the SL(p) magnet with the Lax operator with only one 
diagonal entry containing A. From the results of s.4.3 we can conclude that it has the form 



-(1) 



j{12) 



-i'pxp(A) 



(12) 

-I- 

(13) 



-(13) 



A 



;,(l,p-l) 



V 5. 



(i,p-i) 



(5.11) 



'^^Note that, following [Q, one should consider 2n x 2n Lax matrix to include the UV-finite case (i.e. to 
obtain Q(A) of degree 2n) into the play. Thus doing, one is led to consider also any pn x pn Lax matrices - 
the result we can not interpret. 
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where A is an ordinary constant p x p matrix of rank p — 2 (thus, with zero determinant), detLpxp = 1 and 
the Poisson brackets are 



This Lax operator leads to the spectral curve 



(1,0) 



(1) 



+ P^'>{\)w''-' + ...+Pi''\\)w 



+ 



(5.12) 



(5.13) 



and, therefore, has no n x n representation (in accordance with the argument in the beginning of this subsec- 



tion). However, one can restrict the matrix A further so that the polynomial P^''\\) in (5.13) would be of 
degree n — 2k. This system already has an n x n representation. It can be reached, for instance, in the case 
of p = 3 by the choice 



A = 



1 -1 



in the case of p = 4 by the choice 



/I -1 V 
A^il -1 1 
\1 -1 0, 

etc. In the Appendix we consider how it happens in detail in the p = 3 case. 

The Lax operator in this n x n representation has the following structure (for n > p): 



(5.14) 



(5.15) 



\ 





w 



c-(l) 
^0,p-l 



0,p 



(5.16) 



so that the p — 1 right corner diagonals are filled up by non-zero entries proportional to — and p — 1 diagonals 



below the main one are also non-zero. In formula (5.16) we use the asterisque to denote the non-zero entries 
which we do not write down manifestly for the sake of simplicity. 



Here we use the normalization that slightly differs from (5.3). This Lax operator, after inserting into 
(5.4), indeed, leads to the spectral curve (n > p): 



+ 1 = 



(5.17) 



+ 1 = 



(5.18) 



Moreover, for n = p, one even gets the curve 

+ pP(A)wP-i + Pf_\(A)w''-' - 

that is the maximal possible curve following from (^.4|), at given p and n. 

If now one brings all but the very left two NS-branes to infinity, the system looks as an SU{n) system 



with matter hypermultiplets realized via the semi-infinite branes, the Lax operator (5.16) results into that of 



| |38[ |, and the spectral curve takes the form (5.8). Thus, the construction of jSg], indeed, corresponds to the 
reahzation of the matter hypermultiplets via the semi-infinite branes. 



5.4 Two Lax representations: general comments 

Note that the two "dual" Lax representations has been presented in literature in more invariant terms 
in the Gaudin limit of the Toda chain system. In the 2x2 representation, this limit corresponds to the 
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replacement 7^2x2 — > 1 + eLa with small e and leads to the linear Poisson brackets, instead of (^.(j). Following 
jl^, one can associate this replacement with transition to the continuum limit. In terms of the n x n Lax 
matrices, this limit corresponds to the rational limit of the trigonometric Toda r-matrix. 

The relevant description of the Gaudin Lax representations can be presented via the embedding of this 
finite dimensional integrable system into rational coadjoint orbits of loop algebras. It is based on a family 
of the momentum maps from the space M = {F, G £ M"'^'"} of pairs of n x r rectangular matrices, with a 
natural symplectic structure 

uj = tr(dF A dG^) (5.19) 

to the dual of the positive half of the loop algebra gl{r). This algebra is defined as the semi-infinite formal 
loop algebra on gl{r) that consists of elements X{X) — ^i^'' with Xi G gl{r). Algebra gl{r) as the 

vector space has a natural decomposition 

gl{^) = .g7(^)^ e gli^f (5.20) 

into the spaces of r x r matrix polynomials and strictly negative formal power seria in gl{r) respectively. We 
define the pairing between these two factors 

(X(A), F(A)> = Tr {X{\)Y{\))_^ (5.21) 

where the subscript "-1" refers to the coefficient of A~^. 

Now we fix an n X n matrix A whose spectrum is completely inside a disk D and define the group GL{r) 
of GI/(r)-valued functions of A that are holomorphic inside D. 

Each this matrix A defines a symplectic action of GL{r) on M: 

g{X): {F,G)~^{F„G,), g{\) e Gl(t)^ 

{A - Fg-\\) = {A- \)-^ F, + Fhoi, (5-22) 
g{X)G^ {A - A)-i = (A - A)-^ + gL; 
where Fhoi, Gj^^i are holomorphic in D. 



One can easily check using (5.21) |40| that this symplectic action is Hamiltonian, with moment map 
Ja: M ^ [gUr)"^ defined by 

Ja{F,G) = -G"^ {A-\y^ F (5.23) 

Absolutely analogously to this scheme, one can consider the loop algebra gl{r) and its group GL{r) 
and introduce a fixed r x r matrix Y € gl{r) that define the action 

h{w) : [F, G) -» {Fh, Gh) , h{w) e GUn)^ 
h{w)F {Y ~ w)-' ^ Fh {¥ - w)-^ + Fhoi, (5.24) 

(y - w)-^ G^h-^ = (y - w)-^ Gl + gL 

where -Fhoi, G^oi a^re holomorphic inside disk D containing the spectrum Y . As above, this action is Hamil- 
tonian, with moment map Jy ■ M ^ (^l{n) ^ defined by 

Jy{F,G)^F{Y -wy^G"^ (5.25) 

Now let us define two matrices 

L(A) = Y + G'^ [A- A)"' F, 

(5.26) 

C{w) = A + F{Y-wy^G'^ 
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Using the Adler-Kostant-Symes theorem, one now proves that these two matrices are Lax matrices . Each 
of them is the Lax matrix of the classical Gaudin system. They can be also considered as solutions to the 
moment map equation of the Hitchin system on the sphere with n or, respectively, r marked points so that Y 
and A define the monodromy of infinity. 

The equivalence of these two Lax representations is defined by the fact that they define bi-rationally 
equivalent spectral curves. This follows from the identity 

det [A - A) det {Y + (A - A)"^ F - w) = det (Y - w) det [A + F{Y ~ w)'^ - \) (5.27) 

Note that the construction for the Toda case looks more involved, since, in the 2x2 formalism, it is 
described by the quadratic Poisson brackets. This results in the product of the Lax operators along the chain 
(polynomial in A of degree n). Actually we have to perform the transition from the algebra to the group 
level. Put it difl^erently, this is the product over marked points, instead of sum that could be reproduced 
by the larger matrix still linear in the spectral parameter. Analogously, n x n Toda Lax operator has linear 
Poisson brackets with trigonometric r-matrix. This can be again reproduced by the product of two matrices 
(leading to a polynomial quadratic in w) instead of sum over the two marked points. Unfortunately, there 
is no technique developed for the products of matrices. It looks slightly similar to the averaging procedure 
proposed in |^], which is immediate for the linear Poisson brackets and is quite involved for the quadratic 
ones. 

Let us also note that the natural analog of matrices F and G in the Toda case would be the n x 2 matrix 
with the two columns made of {'0n} and {Xn}- 

Let us still say some words of the algebraic structures in the Gaudin system obtained from the Toda chain 
which corresponds to the above construction at r = 2. In the nxn representation, the Lax operator is defined 
as the solution of the moment map equation for the holomorphic SU{n) connection on the sphere with two 
marked points. The Lax operator in this representation contains momentum variables as the diagonal entries. 
This can be interpreted in the framework of the consideration above as the monodromy at infinity. Turn now 
to the 2x2 representation. Its structure clearly indicates that now momenta play the role of the positions of 
the marked points on the punctured sphere with the holomorphic coordinate A that was the eigenvalue of the 
Lax operator in the first representation. The picture can be clearly generalized to the Gaudin system with 
arbitrary number of marked points in the first representation. 

Returning to the brane interpretation behind the dual Lax representations, it can be viewed as follows. 
In the Type IIA, in one representation we "look from" the diagonal of the brane picture which gives rise to 
the 2x2 representation. It is this look that admits the interpretation in terms of an "effective n-monopole" . 
In the other, nxn representation we take another look and obtain SU{n) gauge group from the parallel 
D4-branes. The DO-branes provide the structure with a few monopoles in the SU{n) theory leading to the 
nxn representation. In the IIB picture, we have f/(l) gauge field on NS5 brane and a possible interpretation 
of 2 x 2 representation is discussed in |^^. The other, nxn interpretation, in this case, comes from the 
wrapped branes around the compact surface. 



6 Monopole ^ spin chain correspondence 

It was already noted by different authors [^ , [isl Q that monopole moduli spaces play an important 
role in the 4d SUSY gauge theories. The evident reason for the appearance of monopoles in this setup is the 
hyper-Kahler structure of the monopole moduli spaces which also can be manifested for the moduli space of 
vacua in SYM theories as well as for the phase space of the integrable systems under consideration. The useful 
approach to the description of the monopole moduli space is to introduce monopole spectral curve and look 
at the moduli space of such curves [^. It was shown [^, Q that the monopole spectral curve is related to 
the spectral curve of the so-called Nahm equations giving an infinite-dimensional counterpart of the ADHM 
construction [^, ^]. It was P.Sutcliffe who first used the Nahm equations in the context of the Seiberg-Witten 
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theories and noted that the solutions of the periodic Toda chain are naturaUy identified with the Cp cyclic 
SU (2) monopoles of charge p. 

It is worth noting that Nahm approach suggests an interpretation of the equations of motion of our in- 
tegrable systems. Indeed, it was mentioned in jisj that the Nahm equations (and, therefore, the equations 
of motion of integrable system) are interpreted as the condition of the BPS saturation which keeps some 
supersymmetries. This means that the demand for the vacuum state in the proper u model to be supersym- 
metric gives rise to the equations of motion. Note also that the fermions in the auxiliary Lax linear problem 
(Baker- Akhiezer function) would be related to the fermion zero modes in a monopole background. 

In this section we demonstrate that such a correspondence between monopoles and integrable chains can 
be pushed further and state that the solutions of the spin chains can exhaust all the solutions to the Nahm 
equations. As a byproduct, we obtain a natural generalization of these equations. Physically the proposed 
correspondence implies a duality of some brane pictures. This can be understood in the framework of the 
Diaconescu construction |43|, who generalized the analogous picture for instantons }47|. 



6.1 S'[/(2)-monopole versus spin chains 

We start with the SU (2) monopoleQ with charge p. It can be described by the Nahm equations that are 
three equations for three p x p matrices Tlit), i = 1, 2, 3 

f (6.1) 

given on the real segment t £ [—1, 1] and regular on the interval t G ( — 1, 1). These matrices satisfy additional 
restrictions 

-7-(t) = -7-t(i), %{t)=T^'{~t) (6.2) 

and boundary conditions 

where 7^^^' ^' form an irreducible p-dimensional representation of SU{2). The numbers 1, -1 are just the 
chosen normalization of the asymptotics of the two Higgs field components at infinity. The physical monopole 
fields are constructed from the zero modes of the linear first order differential operator in the background field 
of the Nahm matrices Ti{t) U, H. 



The three Nahm equations (3.1) can be rewritten as the single equation depending on the spectral pa- 
rameter. For doing this, one needs to introduce the Lax operator 

Ln(\) ^T++To\ + T-X' (6.4) 

where we denoted T± = Ti ± iT2, % = -2iT3. Then, (|t|) can be written as the single Lax equation 

^^ = [A(A),L^(A)], A{X) = l{T+\-'-T-\) (6.5) 

With the Lax representation, one can consider the isospectral problem and construct conservation laws from 
the spectral curve: 

detnxn (Ljv(A)+H=0 (6.6) 

This curve has the form [^, ^ 

w'' + qi{X)w''-' + ... + q^,{X)w''-'' + ... + qp{\)=0 (6.7) 

where (j'fc(A) is a polynomial of degree 2k, and, therefore, the genus of the curve is {p— 1)^ (in accordance with 
(4.8). This curve is an algebraic curve in the space of the oriented geodesies in isomorphic to the complex 



tangent bundle TP^ of the Riemann sphere P^(C) |43|. If x is the local coordinate on this sphere, {w, A) are 



'^Hereafter, speaking of monopole we mean BPS monopole. 
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the standard local coordinates on TP^ given by (w, A) — > w^\x=\. Now conditions ( |6.2| ) and (6^) can be 
rewritten in terms of the spectral curve S ^] as the reality condition 



q^{\) = (-lY\'\,{-\-^) (6.8) 



and the following restrictions: 



is trivial on the curve, L{p — 1) is real (6-9) 

i?°(S,i'"(p- 2)) = for a; G (-1,1) (6.10) 

where L'^ is the linear holomorphic bundle over TP^ given by the transition function exp(— cjw/A) defined on 
the domain A 7^ 0, 00, and L{p) = L (g) 0{p). 

In what follows we establish the "formal" correspondence between Nahm equations and spin chains ignor- 



ing the additional requirements (3.2), (6.3) or (|6.8|)-(6.10) that just restrict the admissible space of solutions 



Thus, we concentrate on the general integrable properties of the Nahm system with no respect to the boundary 
conditions. In fact, these additional requirements are quite restrictive, since the monopole moduli space has 
the dimension 4p increasing with p linearly in contrast to the genus [p — 1)^ and the quadratic grow of the 



dimension of the moduli space of the general curves (6.7) 



First of all, let us comment interpretation of Toda integrable equations as a particular Nahm equation 
(for the symmetric monopole configuration) |42| The key property here is that both coordinates and 

momenta of the integrable systems can be expressed in the monopole terms. More concretely, the Lax matrix 
of the integrable system L(A) can be written as a combination of the Nahm matrices for n symmetric 5(7(2) 
monopoles (see below). However, first we see how the monopoles are getting involved in the picture. Actually 
as we discussed in section 2 we would like to relate the coordinates of Toda particles with the positions of the 
well separated monopoles in a;^-coordinate. The monopoles are localized on D4-branes, one per each brane, 
so we can attribute the same coordinate pi — a;:*i+i5,i34i ^j^^ monopole localized on the i-th D4 brane. Note 
that this is in agreement with the "inversion formula" of the Nahm construction 

T!'"'{t)= /d^x*+'''(a;,t)a;,*'"(a;,t) (6.11) 



which gives the Nahm matrices Ti{t) as the matrices of the coordinates in the Heisenberg picture in the space of 
the fermionic zero modes (fc = 1, .., n) in the monopole background. The relevant Nahm matrices provide 
the momenta pi in the Cartan sector, and coordinates e'^" for positive and negative roots 

where E±j are the generators corresponding to the simple and longest roots of the affine algebra SL{n) and Hj 
are its Cartan generators. Note that the formulae above can be also interpreted in terms of the fundamental 
monopole of SU{n). 

Now let us return to the general monopole (Nahm) system and, first, let us note that the spectral curve 



(6.7) coincides with the spectral curve of the SL{p) spin chain at two sites (4.16). In particular, one can easily 
get from the genus formula ( |4.5| ) for n = 2 the genus (p — 1)^ that is well-known answer for the monopole 
spectral curve Therefore, it is naturally to suggest that these two systems are equivalent. In fact, 

the identity of the spectral curves is still not sufficient to identify two systems. Thus, one needs to check 
the identity of the Lax representations. As soon as there is a gauge transformation that connects two Lax 
operators, the two systems are equivalent possessing the same integrals of motion. In this concrete situation 



should connect the Nahm Lax operator (6.4) and the Lax operator (4.9) of the SL{p) spin chain. In fact, we 
assert that the 2-site transfer matrix 

r2(A) =L2(A)Li(A) (6.13) 



reproduces (6.4). Indeed, the Nahm Lax operator is defined just as an arbitrary matrix polynomial quadratic 



in A, i.e. all the three coefficients {p x p matrices) are arbitrary. On the other hand, the 2-site transfer matrix 
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a) b) 

Figure 5: Two looks at the Toda chain: explanation of the embedding Toda ^ Nahm 



(6.13) is the product of two arbitrary linear matrix polynomials with the unit leading coefhcient. This can 
give us only two arbitrary pxp matrices. The third one comes from the possibihty to multiply spin chain Lax 
operator by any constant matrix without changing the conservation laws (remind that this is the case, since 



the r-matrix ( 4.14 ) is the permutation operator and, therefore, commutes with the expressions like U x U). 



With the established correspondence between SL(p) spin magnet on two sites and the Nahm equations. 



one can immediately explain the embedding (6.12) of the Toda Lax operator into the Nahm one. Indeed, let 
us look at the SU{p) pure gauge supersymmetric theory that is described by the brane configuration of Fig. 5a 
corresponding to the SL{2 Toda chain {SL{2) magnet) given on p sites. This configuration can be equally 
viewed as that depicted in Fig.Sb. This latter point of view implies the interpretation of the same system as 
a specific case of the SL{p) spin magnet given on two sites. This is as stated above equivalent to the Nahm 
system. Since the SL{p) magnet (= Nahm system) corresponding to the brane configuration of Fig.Sb is not 
of the general form, one should expect that the Toda Lax operator (Fig.5a) is only embedded into the Nahm 



Lax operator. This is, indeed, realized in formula (6.12) 



Thus, we have proved that the 5(7(2) p-charged Nahm system is really equivalent to the 2-site SL(p) spin 
chain. This hints that the spin chain at many sites is equivalent to the higher group Nahm systems. In the 
next subsections we check this idea. 



6.2 S'[/(n)-monopole, spin chains and generalized Nahm/Diaconescu con- 
struction 

We start with the description of the SU{p) Nahm construction following [^, In this case, we fix the 
large- distance asymptotics of the n components of the Higgs field to he Hi, i — 1, n, - fii = 0. It naturally 
divides the domain where the Nahm system is given onto n — 1 intervals (/ii, 771^2), (/i„_i,/i„). On each 
a-th interval there is defined the Nahm system for the three pk x Pk matrices Ti^a(t) 

^ = U,,k[T,,a,TkA (6.14) 
Thus, SU{n) monopole is defined by the set of charges (topological numbers) pi, ...,p„_i. 



The system ( 3.14 ) requires some boundary and matching conditions. The details can be found in [ [19| , 
here we only describe the matching condition at the boundary of two intervals /ia for the Nahm matrices with 
some Pa and Pa-i- Let us suppose for definiteness that pa > Pa-i- In the case of equality, all matrices Ti,a{t) 
and Ti,a-i{t) are analytic nearby /la although the condition Ti,a{t) — Ti^a-i{i) around may be not fulfilled. 
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In the case of strict inequality, Ti,a-i{t) is analytic as t ^ fia, while Ti^a{t) has the block form 

%,a-l{ila) + 0{t - ila) 0( - /ia ) ^ ~ ) \ 

J.(Ma) (6.15) 

o{{t-^la)^^--^--^-^^/^) ^^ + o(i) ) 

t- / 

and T^'^"^ realizes an irreducible representation of SU{2). These matching conditions have to be added by the 
boundary conditions po = Pn = 0. 

Now let us turn to the spectral curve. From the described construction it is clear that, by modulo 
boundary and matching conditions, the spectral curve is the product of n — 1 components, each of these being 
the spectral curve for the 5*17(2) monopole of charge Pa- The complete curve, therefore, is a curve of the form 
(6.7) which, in its turn, is a special degenerated curve of the form (4.16) with p = X^^Pa. In the spin chain 
language, this means that we are dealing with a very special configuration of the SL(p) magnet given on the 
n sites. This configuration is determined by "the clasterization" of these n sites into 2-site groups so that there 
is effectively given the SL{pa) magnet on each of these groups, with specific matching conditions. 

To understand the speciality of this configuration better, we need to turn to the Diaconescu construction 
that could allow us to propose a natural generalization of the Nahm construction to the generic SL{p) 
spin chain. 

Looking at this construction, one can easily understand what is the specific of the Nahm system that 
leads to a very peculiar magnet. This specific is obliged to "the locality" of the brane picture in the sense 
that the 1-branes has their endpoints only on the neighbor 3-branes. This implies that the corresponding 
spectral curve is just the product of the curves constructed from each pair of the neighbor 3-branes. In terms 
of spectral curve, this means that there preserved the asymptotic behavior at large w and A, that is w ^ 
Thus, appealing to the spin chain interpretation proposed in this section, one can identify the 3-branes with 
the sites of the spin chain, while the 1-branes gives the SL{p) group, p being the full number of 1-branes. 
This leads us to the natural conjecture that the generic SL{p) spin chain corresponds to the generalized 
Diaconescu picture where the 1-branes connects the very left and the very right 3-branes passing through the 
others (inside) 3-branes. The corresponding Nahm Lax operator is to coincide with the transfer matrix of this 
spin chain (with Lax operators (4.£) twisted by the constant matrices) and can be represented in the form 
ijv(A) = X!r=o'^'^' with n + 1 arbitrary matrices %. 



7 Multiple Aj-scales 

In this section we introduce a new set of gauge models inspired by some generalization of the periodical Toda 
chain which contains many A type parameters. Let us emphasize from the very beginning that the possibility 
to be discussed below is purely nonperturbative. The main idea is to consider the correlated inhomogeneous 
magnetic fluxes and the inhomogeneous "DO brane lattice" that can be obtained by a nontrivial double scaling 
limit from the system described by the elliptic spin Calogero model. 

Thus, we start with the spin generalization of the elliptic Calogero system that can be the most directly 
viewed in the Hamiltonian reduction picture [^, For the simplicity, in this section we discuss the GL{n) 
system, the reduction to the SL(n) case being trivial and immediate. The system is given by the Lax operator 



a{xi — Xj)a{\) 

with the additional requirement Ja — O.Q Thus, we put Jij = (1 — 5ij)fij, fa ~ const. This latter condition 
is necessary for integrability of the system (see fcli). 



'^^Let us note that generally these diagonal elements are just conserved quantities. The choice made is 
consistent with the spin Calogero system obtained from the dynamics of poles of the matrix KP system [pO|. 
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The additional spin variables in the Lax operator (7J) fij satisfy the Poisson brackets 



{fijyfkl} = Sjkfil ~ Silfkj (7.2) 

This symplectic form is generally degenerate and its restriction to the symplectic leaves is related to fixing 
the rank of the matrix fij to be Z < n so that the matrix can be presented in the form 



fa = ^ Uiv] (7.3) 

a=l 



and corresponds to the generic complex orbit GL{n, ([]) standard spinless Calogero system 

GL{n-l,(C) X ((D*)' 

GL(n (T) 

corresponds to the orbit describing the fundamental representation „ , , ' . , , i.e. I — 1. 

Gi(n-1,(D) X (C ) 

The corresponding spectral curve of the spin Calogero system, which generates the conserved quantities, 
is given, as usual |^^, by the determinant form and its genus calculated, say, via Riemann-Roch theorem is 
g = np ~ ^^^^ ' + 1. 

We are not going to discuss here the elliptic system more but instead consider the double scaling procedure 
that provides us with a set of scales in the low-energy theory after its degeneration to the rational limit. Let us 
remind how it works in the spinless case. The steps are as follows jssj. First, we degenerate the bare spectral 
torus T — > ioo. Then, to provide the nearest neighbor interaction, we introduce the homogeneous coordinate 
"lattice" with the large distance A between sites 

Xj = jA + 0, (7.4) 

To see explicitely what kind of interaction emerges in this limit, it suffices to look at the Weirstrass function 
giving the Calogero potential, although the procedure can be easily repeated for the Lax operator. The 
Calogero potential has form 

V{xij) = / ^ p{xij) (7 5) 

i-j 

where Xij = Xi — Xj, = ftjfji is the (coupling) constant that can be made non-depending on indices by the 
proper gauge transformation of fij and the Weierstrass function is defined as 

p(x.,) ^ a^loga(x.,) = ^ -^^^^——^ (7.6) 

m 

In the limit under consideration we introduce the renormalized coupling constant g — goexp{A) so that go 
will be ultimately nothing but Aqcd- Now, choosing A ~ r and taking the limit r — > ioo, we see that only 
m = term survives in the sum (7.6) in the potential (7.5) so that the resulting potential reads as 



Vb(<^0=5oVe^'+^-^' (7.7) 



and describes open (non-periodic) Toda chain. 



In order to get the periodic Toda chain, one needs to fix A = ^. In this case, the m — —1 term in (7.6) 
also contributes into the sum and one finally obtains the potential 

VTci<j}^) = Vo+goe^'-^" (7.8) 

describing the periodic Toda chain. 

Let us now consider the generic spin Calogero system with I = n. In this case, there is a rich spectrum 
of possibilities so that one would expect n scale parameters instead of a single one. Indeed, it is possible now 
to introduce the inhomogeneous "lattice" with n independent parameters which give rise to n infrared scales 
after the independent dimensional transmutation procedure in each (7(1) factor. Quantitatively it corresponds 
to the definition 

Xi+i,i = Ai + (7.9) 
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with arbitrary Ai which are assumed to be large. The potential in this case has the form 

= ^fijfiipi^ij) (7,10) 

Consider now different possible limits. In the simplest case we choose 

fi+i,i = fi,i+i = go,i exp(A,) (7.11) 

with finite go,i assuming that no other fi,j have the same exponential factor. This amounts to the open Toda 
chain with different coupling constants Ai = go^i at each site. Again, in order to get the periodicity of the 
chain, one more constraint has to be imposed 

n 

r = ^Ai (7.f2) 

1=1 

where /„,i = = gio,nexp(A„). This leads to the potential 

n-l 

V{<t^i,) = J2 ahe^'*''^' + 5o,ue^^-^" (7.f 3) 

However, in this simplest case, we get the same periodic Toda chain, since the difference of the coupling 
constants in different sites can be removed by shifts of the Toda coordinates. 

However, there are less trivial limits given by the choice 

fi+3,i = go,i exp(Aj) (7.f4) 



that is to some extent complimentary to (7. If). This choice, with specially adjusted Ai's and the proper 



periodicity condition, leads to a "more filled" Lax operator than Toda one, with several non-zero diagonals 



so that more pairs (j>i,(f)k interact. In fact, this Lax operator is similar to (5.16) and, in the most general 
situation, is "completely filled". This generic system, certainly, gives rise to the system that differ from the 
Toda one but the number of free parameters still can not exceed n (n — 1 in the SL{n) case). 

Turn now to the brane picture behind the described limits. Let us begin with the rough "perturbative" 
configuration. Remind that for the single scale N = 2 theory one has a pair of NS5 brane with n D4 branes 
between them. The distance in the direction is identified with -i^. Therefore, in order to reproduce the 
above consideration, one would order the DO-branes along the a::®-coordinate in the simple Toda case but in 
the general case the DO-branes are to be assigned with arbitrary positions along this direction. This provides 
a set of new scales in field theory. The corresponding multiscale low-energy theory is broken, apart from the 
usual Higgs mechanism of the gauge SU{n) symmetry breaking, by the explicit gauge symmetry breaking at 
the regulator scale. Due to the asymptotical freedom of the theory and the conformal anomaly, this ultraviolet 
breaking is lifted down to the low energy sector. In the Lagrangian terms this means that we break A'' = 4 
theory down to A'^ = 2 adding more complicated regulators with the set of mass parameters that manifestly 
break the gauge symmetry Tr [M($i + "l"!)] > where M is a mass regulator matrix. 

Another interpretation of the regulator scales is related to the recent proposal that the UV cut- 
offs can be identified with the momenta pn in M theory. Therefore, we can describe our multiscale case 
as corresponding to the nondegenerated set of momenta along the eleventh dimension. This implies the 
possibility of nontrivial brane scattering processes with the pn transfer which should has much to do with 
the interpretation of Whitham dynamics as a scattering process [Q . 

Certainly, entering the several UV scales into the low-energy physics deserves further studies and we are 
going to return to this question in the forthcoming paper |5q] . 



8 Discussion on symmetries 

This section is slightly outside the main line of the paper, but we consider it important for the future devel- 
opment. In fact, we speculate on what might be the role of some key notions/properties of integrable systems 
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that have not been used yet in the SUSY theories/brane framework. In particular, we briefly discuss the 
hidden symmetries which are well known as "non-abelian symmetries" in the context of the integrable field 
theories and have some counterpart in the finite dimensional case. Their role in the classical dynamics on the 
spectral curve, is still to be understood but they certainly are of the great importance at the quantum levelQ. 
However, we start with some discussion of how the finite-dimensional system is embedded into 2d one and 
only then come to the symmetry issues. 

8.1 Embedding in 2d integrable systems 

In order to see the full (typically infinite dimensional) symmetry group underlying a finite-dimensional system, 
one needs to embed it into a {2d) infinite dimensional integrable theory. Embedding of such a type is known, 
in particular, for the systems of the Calogero type. In fact this is a heuristic observation that there is a close 
relation of these systems with special solutions in 2d integrable field theories [^. This correspondence has 
no rational explanation yet. 

We expect that the embedding of the many-body systems related to the SUSY/brane theories into 2d 
systems would be very useful in understanding of the processes of creation/annihilation etc. To be precise this 
embedding is formulated as follows . One should start with the rational, trigonometric or elliptic solutions 
to the KP hierarchy with the r-function that has a multiplicative form 



Now the condition for this r-function to satisfy the equation of the KP hierarchy defines the dependence 
of zeros Xi on all the times of the hierarchy. This dependence is determined by the dynamics given by the 
Hamiltonians of the many-body Calogero system, that is to say, by the traces of the Calogero Lax operator 
corresponding to the, accordingly, rational, cylinder or elliptic geometry. The number of factors in the r- 
function coincides with the number of particles in the Calogero system. Two more requirements are to be 
added - the Calogero coupling constant is to be fixed and is not allowed to take arbitrary values and the 
evolution has to be started with the specific "locus" particle configuration with zero particle momenta. This 
correspondence holds for the generalization of the KP system to the 2d Toda lattice hierarchy, the dynamics 
of the specific solutions to the 2d Toda being governed by Hamiltonians of the Ruijsenaars finite- dimensional 
integrable system [ p2[ . 

Let us remark that according to ^] one can establish an analogous correspondence between the Calogero 
and KdV dynamics. In this case, the whole procedure can be naturally performed in two steps. First, one fix 
the system of n coinciding particles with the total charge n, which further decays (in accordance with the KdV 
dynamics) into the zero momenta "locus" configuration that evolves with the Calogero Hamiltonians. This 
locus configuration is not arbitrary but possesses 2n"Sy™™s^ry. This procedure can be also interpreted in 
brane terms, with the Calogero degrees of freedom regarded as DO coordinates. Then, the initial configuration 
is interpreted as a bound state of n DO-branes which KdV-evolve to the stable "locus" configuration which 
then evolves with the Calogero Hamiltonians. 

Within the described picture, the possible role that 2d integrable systems like KP or Toda could play 
looks very promising. Indeed, we have seen that the dynamics of the finite number of DO-branes corresponds 
to the evolution of the very special solution of 2d theory. Therefore, it seems unavoidable that the general 2d 
dynamics is described by the nontrivial DO distribution along the direction. Thus, the integrable system to 
some extend looks as a second quantized "brane or string" theory. Note that the KdV dynamical variable u(x) 
now can be treated as the a;''-component of the stress energy tensor so that the expression for the simplest 
(rational) solution 



Let us fix at the very beginning that all the hidden symmetries under consideration are the Poisson-Lie 
ones. 



r{x,t) = Y[ (T{x-Xi{t)) 



(8.1) 
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in DO brane terms looks quite satisfactory. 

This second quantized picture can be pushed even further. To this end, let us note that embedding 
different Calogero systems into the same KdV (or, more generally, KP/2d Toda lattice) hierarchy plays a 
unifying role. Indeed, solutions of the Calogero systems with different number of particles, that is to say, 
with different gauge group ranks, being considered as solutions to the same KdV hierarchy are related by a 
Backlund transformation that can be realized by the (fermionic) vertex operators |^]. Thus, these vertex 
operators are associated with the operator of the brane creation providing the change of the gauge group 
rank. This means that the (sub)space of the solutions to the larger hierarchy can be chosen as a configuration 
space of all the vacua/string theories. In its particular realization, this space can be described as an infinite 
dimensional Grassmannian ^] and reminds of the analogous construction known from the perturbative string 
physics 

In the rest part of this section we discuss the issues of symmetries in integrable systems that can be easier 
and more naturally formulated for the quantum systems. Although so far we have met only the classical 
integrable systems, below we follow the general approach equally applicable to both quantum and classical 
integrable systems. Among other advantages, this will make quantization of the construction an immediate 
thing. 

Let us give some general comments concerning the quantum picture. From the consideration above it is 
clear that the integrable dynamics would proceed on the instantonic moduli space so all the ingredients of the 
quantum picture should acquire the meaning of the appropriate objects on this moduli space. Note that in 
the course of quantization, i.e. constructing the wave function, in the present "Hitchin treatment" we need to 
neglect the Whitham variables taking them as just parameters. This corresponds to the Born-Oppenheimer 
approximation. It is just in this approximation the whole wave function reflects the structure of the instanton 
moduli space. 

As we have argued, instead of the infinite dimensional instanton moduli space, we work with the finite 
number of effective variables which substitute the infinite instanton sums. The wave functions seem to be 
related to the cohomology of the moduli space while the quantum partition function would play the role of 
the generation function for the intersection numbers in cohomologies. 

Asymptotics of the wave function gives the S-matrix for the effective degrees of freedom - DO branes. 
The integrability and the presence of the additional conservation laws promise the exactness of the S'-matrix 
as well as its pure pole structure. At the same time, in order to derive the form-factors (in our context, it 
corresponds to the form-factor of the compound state of multiple DO branes), one needs to use the hidden 
non-abelian symmetries which we are going to discuss now. 



8.2 Backlund transformation 

The remnant of the infinite dimensional non-abelian symmetry in the Toda system is the Backlund transfor- 
mation. It is known that the nontrivial Backlund transformation for the periodic Toda chain corresponds to 
the transition between two equivalent sets of variables X2n and X2n+i- At the quantum level, it is convenient 
to consider the wave function in the action-angle representation which is the solution to the Baxter equation 

Q{\ + A^gSoQi^ - «) = Q{X)TrT{X) (8.2) 

where T{X) is the transfer matrix of the Toda chain and A is the spectral parameter. This equation actually is 
the quantum counterpart of the spectral curve equation considered as the operator acting on the wave function 
in the separated variables [ p7[p°| . 

The key point is that the function Q turns out to be the generating function for the canonical Backlund 
transformation. Namely, it can be related |5q] to the kernel of the operator acting between two set of variables 



^°In the classical limit, the functions ^^J^^^ and '^q(^x)^ can be identified with w and respectively. 
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qijQi- Explicit form of the integral kernel of Q{X) looks as follows 

j 

where 

W\{q) — exp{iXq — e'') 
_ (8.4) 

W\{q) = exp{—iXq — e') 

Logarithms of the matrix elements of Q are the generating functions for the canonical Backlund transformation. 
This is in agreement with the general viewpoint: Q is actually related to the Toda r-function while the r- 
function is the generating function for some canonical transformations. 



8.3 Semiclassical Yangian type symmetry 

Let us turn now to the analogous symmetry in the integrable system corresponding to the N = 2 SQCD - 
spin chain [^. It is known that the spin chains enjoy two different sets of the conservation laws - Abelian 
and non- Abelian ones (they are local and non-local correspondingly). The charges of the both sets are 
conserved but, because of the non-abelian nature of the second set, one can not diagonalize the Yangian 
charges simultaneously. The explicit formulae for the Yangian structure can be found in [|9[. Hereafter, we 
consider the SL{2) spin chain, it is homogeneous and the Casimirs at all sites are the same. This is equivalent 
to the coincidence of all masses of the fundamental matter in the N = 2 SQCD. 

We now emphasize that the knowledge of charges is equivalent to the data of the monodromy matrix T 
(we use here the unit normalization of the determinant): 

Yangian charges can be expressed in terms of the monodromy matrix as follows 

oo 

8.6 

oo ^ ' 

'P^{X) = ^A-"-Vr = e,,kVjk{X) 1,2,3 

n = 

Therefore, the generators Vij{X) and Vi{X) are quadratic functions of the matrix elements of T{X). 
It is also possible to express r(A) in terms of the generating functions Vi{X): 

T(A) = \W(X)1 - '-W-\X) J2 ■P^{X)a, (8.7) 



with W{X) = \/2 + ^/4~P^{X). 



Thus, the charges Pi (A) contain the same amount of information as the monodromy matrix. There are 
also the following relations between the "Pi (A) and the matrix elements of r(A): 

P+(A) = Pi(A) + ip2(A) = 2iW(X) C(A) 

V-{X) = Pi (A) - ip2(A) = 2iW(X) B(X) (8.8) 
Trr(A) = W{X) = W [P^(A)] 
Therefore, P^(A) is also a generating function for (non-local) commuting quantities. 
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For the completeness, let us manifestly write down the Poisson brackets of semiclassical Yangian generators 

(8.9) 

{{vlvj}, {vl Pl}} + {{vlPl}, {V?, p}}} = 8(A™76fe,„ + 

Wltn /ijjj. — g^ila^jmb^knc^ 

One sees that these relations form a deformation of those defining the Borel part of the su{2) loop algebra. 
All the non-local charges P"j can be expressed in terms of the Poisson brackets between the two first charges 
Pi and Pl- Therefore, the whole Poisson algebra of the symmetries is generated by these two charges. They 
are nothing but the corresponding Chevalley generators. 

Finally, let us describe how the monodromy matrix generates the transformations of the spin variable 
on the fc-th site which corresponds in the brane language to the action on the fc-th (D6+D4) brane. By an 
exphcit computation of the Poisson brackets between the monodromy matrix and the spin variables, one can 
show that the variation S"S{k) of the spin variables reads as: 

5:Sik) = {pr,Sik)}^ */|^^" ^^^{ ('^»r-'(A)® l){r(A)®l,l®5(fc)} ) (8.10) 

Here Tri denotes the trace over the first space in the tensor product and Vi is the parameter of the transfor- 
mation. This formula describes the transformations S{k) — > 5"S{k) and indicates that r(A) also generates 
the non-Abelian symmetry. 

Let us comment now the possible meaning of the non-abelian symmetry within the brane approach. 
Usually this symmetry is important at the quantum level providing the set of restrictions onto the S'-matrix. 
In our context it would mean that the symmetry restricts the brane scattering. Note that the Yangian charges 
are essentially non-local. In the SQCD context this would correspond to mixing the different flavors (that 
are associated with different sites of chain) and presumably might be relevant for the barionic branch of the 
theory. 



9 Conclusion 

In the present paper we discussed effective N = 2 SUSY field theories, corresponding integrable systems, 
monopoles and their stringy analogues. The equivalences were established at low energies in field theories 
as well as in string theory models. One can try to go further and ask whether such connections persist if 
we consider higher derivatives in field theory corresponding to non-lowest excitations in string theory. The 
answer seems to be yes, and there is natural candidates in the world of integrable systems for the counterpart 
of higher-derivative terms. They are related to the higher times and introducing two-dimensional integrable 
systems. We are going to elaborate this point in the subsequent publication jsGj. We also postponed the 
discussion of elliptic models (and their degenerations) as well as five- and six-dimensional field theories (XXZ 
and XYZ spin chains), which constitute the two last rows in Fig.l. The new examples confirm the validity 
of the general approach and suggest the proper way to include additional degrees of freedom. Moreover, 
it provides the possibility to get completely new viewpoints in the field theory setup which we have seen 
in section 7. Nevertheless, despite many supporting arguments including the identification of branes as the 
proper degrees of freedom, the clear physical answer to the question: "Why the low-energy SYM is governed 
by integrable system?" is still missed. 

The answer to this question requires the knowledge of the exact derivation of the effective degrees of 
freedom in the integrable systems from the instantonic sums. We have argued that somehow the instantons 
are summed into the finite number of degrees of freedom. However, the proper derivation and explanation of 
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this phenomena from the path integral of SYM is still lacking. The answer would imply a nontrivial localization 
while integrating over the instanton moduli space. 

So far only the tiny part of the "integrability world" has been recognized in the supersymmetric YM 
theories. Say, the issues of underlying symmetries especially in connection with the quantization problem as 
well as interpretation of the non-Abelian "spectrum generating" symmetries in the brane language have been 
out of mainstream. We are planning to discuss this points elsewhere. 

Another important point about the integrable treatment is to interpret the wave functions (which would 
be interpreted as the correlators in a topological theory) that should provide the important information 
about the vacuum configurations. The wave functions actually serve as generating functions for the canonical 
transformations. It is interesting to identify these canonical transformations in the brane language. 

Let us say some words about the possible role of integrability in TV = 1 SUSY theory. It was recently 
shown [ ^ , |6^ ] that the perturbative spectral curve in such a theory is a sphere with marked points which 
suggests that we leave the family of the "finite-gap" solutions. But it can be not the whole story. Indeed, 
the branching points on the spectral curves come from the Coulomb branch moduli and so Aqcd does. The 
Coulomb moduli degenerate to the set of points which means that from the N = 2 SUSY point of view we 
can no longer vary the integrals of motion. But the potential cut in the strong coupling region still persists 
due to Aqcd so the one-gap solutions are potentially possible. 

Recently another object from the integrability cuisine was also discussed in the context of = 1 theory, 
that is, the "creek equation" from j67j that gives the profile of the BPS objects. This equation is a direct 
counterpart of the Nahm equation treatment above. So the would be integrable system in A" = 1 context 
corresponds to the dynamics of several interacting BPS states . 

We have seen some similarities with the "confining string" scenario above. Indeed, one of the main points 
of that scenario is the presence of "the monopole ring" in the vacuum state, the appearance of the effective 
"confining string" to substitute the infinite number of instanton and the appearance of some auxiliary surface 
whose genus corresponds to the rank of the gauge group. One immediately recognizes these ingredients in the 
integrability approach but we do not know any real statements behind this observation. 

Note that the picture considered in the paper has, in fact, some solid state analogies. There are many 
examples in the solid state physics when the mass gap in the spectrum of the quasiparticles vanishes on (hy- 
per)surfaces of different dimensions in momentum space. Just these "defects" ( "monopoles" ) in the momentum 
space are counterparts of branes in our approach. The very appearance of the (topologically nontrivial) hy- 
persurfaces in the theories with anomalies has a simple interpretation. Indeed, the presence of the anomaly 
(for instance, the conformal anomaly in the SYM theory) implies the nontrivial level crossing that can be 
seen within the standard Berry phase treatment. In our consideration the mass of regulator (entering the 
low-energy sector due to the conformal anomaly) gives rise to the "monopole field" in the space of fields in 
agreement with this interpretation. 
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10 Appendix. SL{3) chain 

In this Appendix we consider the simple case of the SL{3) chain to give exphcit examples of some formulas of 



>S'L(3) spin chain. The Lax operator is given by the formula 



Li 



s'. 



(1) 



V 



s_ 



(12) 



The spectral curve in this case has the form 

^3 ^ 'prT(A)m^ + ^ detMi(A)M + detr(A) = 



(A.l) 



(A.2) 



where T{X) is, as before, the transfer matrix and Mi{X) is matrix obtained from the transfer matrix by removing 
the i-th column and i-th row. In this reference integrable system, Ji = 1 and J2 = detr(A) = J^^detLi(A). 
To determine mass spectrum we need to calculate detLi(A). It has the form 



I 0(1) 0(2) 



detL,(A) = (A + X^f - C^''(A + A,) - Q'' 

where Casimir functions are 

-^(i) _ c(12) c.(12) c.(2) , c(2) c(2) c.(l) 0(1) c{2) c{12) c{2) c,(12) 

"-'a ^ '-' + ,1 '-'0,i ' 1*^0,1 + '-' - .i'-' ~ .i'-' + ,i ' 

,qW 0(2) f qW , q(2)\ ^(1) ^(1) / ^,(1) c,(2)\ 

that can be manifestly checked using the Poisson brackets 

rqW qW\_q(2) qW f q(2) c{2) 1 _ ^(l) r, q(2) r qW \ _ qW 

1'-'+ I"-"- / — '-'o ~ ''o ' 1''+ I"-"- /— 'Jq ^''o ' '^0 i — ^'^± I 



(A.3) 



(A.4) 



(2) c.(l)l 



^(2) c.(2)l 



;.(2) 



and defining relations 
so that 



;.(2) 

^(12) 



^^(2) ^ ^(12) ^ ^ ^(1) ^ ^^(1) ^ ^(12) ^ ^ ^(2) ^ ^^(2) ^ ^(12) } 



(12) ^(2)-, 



j(12) _g.(12) 



(A.5) 



(A.6) 



(A.7) 



Generally, the determinant (A.3) has three distinct roots (masses), since we deal with the reference system, 
that is to say, Ji(A) = 1. In order to obtain non-unit Ji(A), one needs to consider the case of, two coinciding 
roots. This condition implies specialization of the general orbit 



and the coinciding masses are equal to 



4C| = 27C| 



^ ^ /^2 

mc = -Xi ± \ l — 



(A.8) 



(A.9) 



Then, the factor (A — nic)'^ is nothing but Ji(A) (Ji(A) of the general form can be obtained by degenerating 
at a set of sites). Now one can check that condition (A.5) leads to factoring out the multiplier (A — rric) in 
the coefficient in front of w in (A.S) to reproduce the correct dependence of the spectral curve on Ji(A). This 
coefficient is equal to 

T11T22 + T11T33 + T22T33 — T12T21 — T13T31 — T23T32 (A. 10) 
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For the definiteness, let us consider the degeneration at the first site. Note that ( |A.10|) can be rewritten as 



(I/l,llLl,22 — il,12il,2l) ^riir22 — Ti2T21^ + 

+ (i/i,iiLi,33 — -£/i, 13^1,31) ^riir33 — ri3r3i^ + 

+ (il, 22^/1, 33 — il, 23^1, 32) I r22r33 — 723732 1 



(A.ll) 



where T is the transfer matrix without the first site and -£/i,ij is (i, j)-matrix element of Li. This expression 
can be rewritten as follows 



X + X^ + 



X A + Ai + 



X T11T22 - T12T21 + 



(A.12) 



X Tiir33 - ri3T3l + 



+ (A + A. + is« + i^(^ 



X ^r22r33 — r23r32^ 

Now by the straightforward calculation one can check that each of three terms in this formula is proportional 



to (A — rric) provided the condition {\-.S) is satisfied. 



S'L(3)-Toda-like system. Now let us consider the "limiting" degenerations of the Lax operator (A.l) 
analogous to the Toda system in the SL{2) case. We follow SL{2) case and multiply the Lax operator (A.l) 
by the constant diagonal matrix U . In this case, there are two possibilities. First possibility is to choose 



U = 



/ 1 
a 

V 1 y 



(A.13) 



and redefine S^^^ ■■ ^^-^ ^ i^^-^ ■ Now, bringing a to zero, one obtains the Lax operator (hereafter, 

we omit inhomogeneity, since it does not eflect the result) 

^(12) ^ 





( 5^^^ + A 


5(1) 


L = 














j(2) 



(A. 14) 



+A J 
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Then, instead of (A^), we get 



and instead of (A.C) 



(A.15) 



(A.16) 



while (A.7) remain unchanged (now, however, one can not use the defining relations (A. 6) to generate all other 
(A. 7) but they should be given independently). 



At the next step, the second Casimir function of the algebra ( A.15 )-( A.16 ) S_ + S_^_ S_ should be 
put equal to zero. Then the determinant of the Lax operator (A. 14) is equal to the third Casimir function 

^(1) 5(2)5(12) ^ 5(1)5(2)5(12) ^ 5(1)5(1) (^2Si'^ + 5(2)^ (A.17) 
that can be chosen equal to unity. 



The second possibility of degenerations of the Lax operator (A.l) looks slightly simpler having less non- 
trivial entries. It can be achieved by multiplying the Lax operator by the matrix U of the form 



(A.18) 





/ 1 







u = 





a 







[ 





a ) 



This form of the matrix allows one to redefine, apart from 5*' 



(1) 



1 q(1) q(2) 



15(2) 5(12) 



^Si''\ also 

a -r ' 



the diagonal element S^'' a^(^^ without tuning the inhomogeneity (compare with the SL{2) case). Now 
taking a to zero, one obtains the Lax operator 



(A.19) 







5(1) 


5(12) X 


L = 




q(2) 


5(2) 











and algebra of the Poisson brackets 



;(2) 



;(2) q{2) 



(1) c.(l)l 



j(l) 



j(l) c{2)i 



(1) c{2)l 



-jW c(12)i 



j(2) 



(1) q(12)| ^ 5(2) {5^2) q(1)i _ 4_e{12) /c'(12) q(2) 



{sx',s: 



(2) 

^5(2)^^(12)^ ^15(2)^ 5(12)} ^0 
(12) c,(12)l _ oqW , q(2) 



(A.20) 



One can easily check that Sq^"* and S^^' are the Casimir functions of the algebra (A.20). In order to get 
constant determinant of the Lax operator (A.19) we put S^^^ = ±1 and 5*0^' — 1. Finally this leads us to the 
Lax operator 



L = 





5(1) 




1 




1 



(12) 



-1 / 



and the Poisson bracket algebra 



The determinant of the Lax operator (A.21) is equal to the Casimir function 

Q< ^ 5(1) 5(1) ^ 5(12) _ 5(^12) 5(12) _ ^(1) ^(12) 

and also can be put equal to unity. 



(A.21) 



(A.22) 



(A.23) 
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Let us note that ( A. 21 ) describes a slightly specified case of the general situation, since the constant 
/I -1\ 

sub-matrix I ^ ) ™ operator is chosen to be traceless. It can be avoided by a proper rescaling 

of S'g^' and inhomogeneity. This particular choice is one that admits the n x n representation. As any n x n 
representation for the SL{p), p > 2 magnet, it describes a degenerated spectral curve as compared to the 
general pure gauge theory (see s.5.3). Let us trace out how it happens in this concrete case. 



First, one can observe that the system of the Poisson brackets ( A.22| ) has Casimir functions =F S'^^\ 



Taking into account that we put Casimir function ( A. 23) to be unit, we can choose Sj^ 



j(12) 



(12) 



T c±, 



c+c_ — 1. The simplest choice is c+ — = 1. Then, literally following s.5.1, one gets the nxn representation: 



'-'0,1 

C{1) , qW 

'-'+,2 + '-'-,2 



•^0,2 

q(1) , 

'-'+,3 + ^-,3 



c.(l) 
'-'0,3 



V 



'-'O.n 



(A.24) 



(1) q(1) 



with the following Poisson brackets: {S^i. ^,5* 
two harmonic oscillators {pa, qb} = Sat 



} 



a, 6= 1,2: 5'''' = ple''^ Si'' = qie 



_i_ jUij that can be realized by the 

-12 



The Lax operator (A.24) leads to the following spectral curve (if n > 2) 



-2(^)- + — 

W 







(A.25) 



In fact, the low diagonal of units in the Lax operator is nothing but the product of c+c_. Requiring this 
product to be zero, we can reproduce the particular Lax operator of the construction (see s.5.2) and the 
spectral curve takes the form 

u; + pW(A) + PfJ2(A)i = (A.26) 
Physically this corresponds to bringing one of the NS5-branes to infinity that results to the two NS5-branes 



remaining with the semi-infinite branes attached to one of these NS5-branes (cf. (A.26) with (5.10)). 

To conclude this Appendix, we see that there are, indeed, two different degenerations, one being charac- 



terized by the Lax operator (A. 14) that contains two diagonal elements with the spectral parameter A, and 



the other one whose Lax operator contains only one diagonal element with A. Let us note that, unlike the 
SL{2) case, bosonization of the both 5*1/(3) limiting degenerations looks quite involved. 
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